 Prefaţă Cartea de faţă fost elaborată în cadrul proiectului POSDRU/ /S/ , “Formarea cadrelor didactice universitare şi studenţilor în domeniul utilizării unor instrumente moderne de predare-învăţare-evaluare pentru disciplinele matematice, în vederea creării de competenţe performante şi practice pentru piaţa muncii” Finanţat din Fondul Social European şi implementat de către Ministerul Educaţiei, Cercetării, Tineretului şi Sportului, în colaborare cu The Red Point, Oameni şi Companii, Universitatea din Bucureşti, Universitatea Tehnică de Construcţii din Bucureşti, Universitatea „Politehnica” din Bucureşti, Universitatea din Piteşti, Universitatea Tehnică „Gheorghe Asachi” din Iaşi, Universitatea de Vest din Timişoara, Universitatea „Dunărea de Jos” din Galaţi, Universitatea Tehnică din Cluj-Napoca, Universitatea Decembrie din Alba-Iulia, proiectul contribuie în mod direct la realizarea obiectivului general al Programului Operaţional Sectorial de Dezvoltare a Resurselor Umane POSDRU şi se înscrie în domeniul major de intervenţie Calitate în învăţământul superior Proiectul are ca obiectiv adaptarea programelor de studii ale disciplinelor matematice la cerinţele pieţei muncii şi crearea de mecanisme şi instrumente de extindere oportunitãţilor de învãţare Evaluarea nevoilor educaţionale obiective ale cadrelor didactice şi studenţilor legate de utilizarea matematicii în învăţământul superior, masterate şi doctorate precum şi analizarea eficacităţii şi relevanţei curriculelor actuale la nivel de performanţă şi eficienţă, în vederea dezvoltării de cunoştinţe şi competenţe pentru studenţii care învaţă discipline matematice în universităţi, reprezintă obiective specifice de interes în cadrul proiectului Dezvoltarea şi armonizarea curriculelor universitare ale disciplinelor matematice, conform exigenţelor de pe piaţa muncii, elaborarea şi implementarea unui program de formare cadrelor didactice şi studenţilor interesaţi din universităţile Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie partenere, bazat pe dezvoltarea şi armonizarea de curriculum, crearea unei baze de resurse inovative, moderne şi funcţionale pentru predarea-învăţarea-evaluarea în disciplinele matematice pentru învăţământul universitar sunt obiectivele specifice care au ca raspuns materialul de faţă Formarea de competenţe cheie de matematică şi informatică presupune crearea de abilităţi de care fiecare individ are nevoie pentru dezvoltarea personală, incluziune socială şi inserţie pe piaţa muncii Se poate constata însă că programele disciplinelor de matematică nu au întotdeauna în vedere identificarea şi sprijinirea elevilor şi studenţilor potenţial talentaţi la matematică Totuşi, studiul matematicii evoluat în exigenţe până a ajunge să accepte provocarea de folosi noile tehnologii în procesul de predareînvăţare- evaluare pentru face matematica mai atractivă În acest context, analiza flexibilităţii curriculei, însoţită de analiza metodelor şi instrumentelor folosite pentru identificarea şi motivarea studenţilor talentaţi la matematică ar putea răspunde deopotrivă cerinţelor de masă, cât şi celor de elită Viziunea pe termen lung acestui proiect preconizează determinarea unor schimbări în abordarea fenomenului matematic pe mai multe planuri: informarea unui număr cât mai mare de membri ai societăţii în legătură cu rolul şi locul matematicii în educaţia de bază în instrucţie şi în descoperirile ştiinţifice menite să îmbunătăţească calitatea vieţii, inclusiv popularizarea unor mari descoperiri tehnice, şi nu numai, în care matematica cea mai avansată jucat un rol hotărâtor De asemenea, se urmăreşte evidenţierea noi motivaţii solide pentru învăţarea şi studiul matematicii la nivelele de bază şi la nivel de performanţă; stimularea creativităţii şi formarea la viitorii cercetători matematicieni unei atitudini deschise faţă de însuşirea aspectelor specifice din alte ştiinţe, în scopul participării cu succes în echipe mixte de cercetare sau abordării unei cercetări inter şi multi disciplinare; identificarea unor forme de pregătire adecvată de matematică pentru viitorii studenţi ai disciplinelor matematice, în scopul utilizării la nivel de performanţă aparatului matematic în construirea unei cariere profesionale Am încercat să facem cât mai atractivă şi accesibilă prezentarea, simplificând expunerea fără pierde din rigoarea matematică rezultatelor Lucrarea este structurată în patru capitole, ultimul referindu-se la probleme de stabilitate clasică şi urmăreşte în principal subiectele prevăzute în programa actuală de studiu, cu precădere cele care pot servi la rezolvarea problemelor tipic inginereşti Astfel, fiecare capitol se încheie cu un paragraf de aplicaţii în diverse domenii: mecanică, astronomie, hidrotehnică, statica construcţiilor, etc Sunt modelate probleme concrete simple, folosind ecuaţii diferenţiale ordinare Prezentarea aplicaţiilor este realizată în patru etape: problemă fizică, model matematic, determinarea soluţiei şi interpretarea ei fizică Considerăm că numeroasele legături cu disciplinele inginereşti, legături pe care le-am pus în evidenţă prin aceste aplicaţii, fac cu atât mai convingător studiul ecuaţiilor diferenţiale ordinare pentru studenţii din universităţile tehnice Paragrafele însoţite cu asterisc pot fi omise, ca şi serie de demonstraţii Le-am introdus, totuşi, pentru unitatea şi logica expunerii Menţionăm că ele sunt, de fapt, destinate studenţilor celor mai interesaţi de domeniul ecuaţiilor diferenţiale şi care văd în viitoarea lor profesiune nu numai un mijloc de trai, dar şi cheie esenţei fenomenelor naturii; ei caută cu perseverenţă “sâmburele” matematic care guvernează din abstract aceste fenomene, căci doar el asigură viziune completă şi unitară asupra fenomenelor studiate şi, deci, prevederea şi stăpânirea acestora Conţinutul teoretic al primelor trei capitole fost realizat de prof Ileana Toma şi conf Emil Popescu, de la Universitatea Tehnică de Construcţii din Bucureşti, iar cel al capitolului – de conf Aurelian Cernea de la Universitatea Bucureşti Aplicaţiile în mecanică şi fizică au fost realizate de conf Dan Comănescu şi conf Ioan Caşu de la Universitatea de Vest din Timişoara, precum şi de echipa Universităţii Politehnice din Cluj, formată din conf Gloria Cosovici şi conf Sorin Comşa Aplicaţiile în mecanica construcţiilor aparţin regretatului profesor M V Soare şi au fost publicate în cadrul Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie volumul “Ecuaţii diferenţiale cu aplicaţii în mecanica construcţiilor”, tradus în Springer coautori: P P Teodorescu, Ileana Toma Bibliografia cuprinde şi link-uri cu site-uri pe care studenţii pot consulta şi online manuale cuprinzând tematici de ecuaţii diferenţiale ordinare Autorii CUPRINS PREFAŢĂ CAPITOLUL ECUAŢII DIFERENŢIALE ORDINARE DE ORDINUL ÎNTÂI Noţiuni preliminare Exemple Formele sub care se prezintă ecuaţiile de ordinul şi soluţiile lor Forme ale ecuaţiilor de ordinul I Forme ale soluţiilor Tipuri de ecuaţii diferenţiale de ordinul rezolvabile prin cuadraturi Ecuaţii cu variabile separate Ecuaţii cu variabile separabile Ecuaţii diferenţiale omogene, de gradul m Ecuaţii cu diferenţiale totale exacte Factor integrant Ecuaţii diferenţiale lineare de ordinul I Ecuaţia Bernoulli Ecuaţia Riccati Ecuaţia Clairaut Ecuaţia Lagrange Metoda aproximaţiilor succesive Teorema clasică de existenţă şi unicitate Cauchy-Picard Principiul contracţiei Aplicaţii în mecanică, fizică şi inginerie CAPITOLUL ECUAŢII DIFERENŢIALE ORDINARE LINEARE, DE ORDINUL Noţiuni preliminare Exemple Ecuaţii diferenţiale lineare şi omogene de ordinul Ecuaţii diferenţiale de ordinul lineare şi neomogene Ecuaţii diferenţiale lineare de ordinul cu coeficienţi constanţi Ecuaţii diferenţiale lineare şi omogene Polinom diferenţial Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Ecuaţii diferenţiale lineare şi neomogene Ecuaţii diferenţiale de ordin superior, integrabile prin cuadraturi Ecuaţii reductibile la EDO cu coeficienţi constanţi Aplicaţii în mecanică, fizică şi inginerie CAPITOLUL SISTEME DE ECUAŢII DIFERENŢIALE ORDINARE Sisteme de EDO de ordinul lineare Sisteme de EDO de ordinul lineare, cu coeficienţi constanţi Exprimarea soluţiei unui sistem de EDO lineare folosind exponenţiala de matrice Sisteme de ordinul nelineare Sisteme simetrice Integrale prime Aplicaţii în mecanică, fizică şi inginerie CAPITOLUL STABILITATE Stabilitatea soluţiilor ecuaţiilor diferenţiale Stabilitatea Liapunov Funcţia Liapunov Sisteme dinamice autonome Comportament pe termen lung al soluţiilor Aplicaţii în mecanică, fizică şi inginerie REFERINŢE BIBLIOGRAFICE CAPITOLUL ECUAŢII DIFERENŢIALE ORDINARE DE ORDINUL ÎNTÂI NOŢIUNI PRELIMINARE EXEMPLE Se ştie ce este aceea ecuaţie algebrică ecuaţie diferenţială este şi ea o egalitate, ce admite însă ca necunoscută funcţie şi mai cuprinde şi derivatele acesteia Deosebim două posibilităţi: aplicaţii  funcţia necunoscută depinde de singură variabilă şi atunci vom avea o ecuaţie diferenţială ordinară prescurtat EDO Aplicaţii  funcţia necunoscută depinde de mai multe variabile, caz în care vom avea o ecuaţie cu derivate parţiale prescurtat EDP Subiectele tratate în cadrul acestui curs aparţin cazului Forma generală unei ecuaţii diferenţiale ordinare este, conform celor spuse anterior, F Definiţia Numim ordin al unei ecuaţii diferenţiale ordinare ordinul maxim de derivare al funcţiei necunoscute y Una dintre problemele esenţiale ale calculului diferenţial este aceea de a determina derivata unei funcţii date Cea mai simplă problemă inversă aparţine calculului integral: Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie PROBLEMĂ Dându-se funcţie reală, de variabilă reală, să se determine primitiva sa Dacă notăm primitiva lui cu atunci formularea matematică acestei probleme este: f x y = d d sau, echivalent dy dx Relaţiile de mai sus sunt, de fapt, cele mai simple ecuaţii diferenţiale şi ştim cum să le rezolvăm Într-adevăr, ştim că cea mai generală funcţie satisfăcând sau este y ∫ dx O primitivă arbitrară lui poate fi deci numită soluţie ecuaţiei Introdusă în ea conduce la identitate Deci şi în cazul ecuaţiilor diferenţiale, soluţie transformă ecuaţia într-o identitate, exact ca în cazul ecuaţiilor algebrice În expresia semnul desemnează una dintre primitivele lui iar este o constantă arbitrară Deci funcţia nu este determinată în mod unic de ecuaţia sau astfel încât putem spune că ele admit infinitate de soluţii Fiecare din aceste soluţii se pot determina dând lui diferite valori numerice Terminologie § Soluţia ecuaţiei se numeşte soluţie generală § Orice soluţie obţinută din soluţia generală prin particularizarea constantei se numeşte soluţie particulară Ecuaţii diferenţiale ordinare de ordinul întâi § soluţie care nu se obţine din cea generală prin particularizarea constantei arbitrare se numeşte soluţie singulară După toate aceste consideraţii, s-ar părea, la prima vedere, că ecuaţiile diferenţiale au apărut într-un cadru strict matematic, ca completare logică formală calculului diferenţial Acest domeniu al matematicii îşi are însă originea istorică în mecanica newtoniană Newton, iniţiatorul calculului diferenţial alături de Leibniz, modelat cu o surprinzătoare intuiţie serie de fenomene fizice prin ecuaţii diferenţiale Astfel, faimoasa lege II-a mecanicii enunţată pe scurt: “Rezultanta forţelor ce acţionează asupra unui sistem este egală cu produsul dintre masa sistemului şi acceleraţia acestuia”, lege care, de altfel îi poartă şi numele, se exprimă matematic sub forma: ma şi nu reprezintă altceva decât un sistem de ecuaţii diferenţiale Într-adevăr, acceleraţia este derivata doua deplasării în raport cu timpul; această observaţie aparţine unui alt titan al ştiinţei, Leonhard Euler Pentru edificare, să urmăm drumul propus de Newton în studiul unui caz foarte simplu PROBLEMĂ Să se studieze mişcarea pe axă verticală unei particule punct material sub acţiunea propriei greutăţi Rezolvare Construim mai întâi modelul matematic Trebuie deci să determinăm a funcţia necunoscută funcţiile necunoscute cărei cunoaştere înseamnă cunoaşterea fenomenului; b legea fizică legile fizice care guvernează fenomenul Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Presupunem că Oy este axa verticală de-a lungul căreia cade particula, originea fiind situată la suprafaţa pământului vezi figura de mai jos Mişcarea particulei este cunoscută dacă se cunoaşte unica sa coordonată anume, poziţia sa pe axa Oy în fiecare moment t Funcţia necunoscută problemei este deci cu semnificaţia fizică de deplasare particulei În problemele de mişcare, legea a doua lui Newton joacă un rol esenţial Aplicând-o pentru unica componentă acceleraţiei, găsim ma mg m fiind masa particulei iar modulul acceleraţiei gravitaţiei Semnul minus provine din faptul că axa Oy este dirijată în sus, iar forţa de gravitaţie în jos Ţinând seama că acceleraţia este derivata doua deplasării în raport cu timpul şi simplificând cu m, rezultă g t y d d Ecuaţia reprezintă modelul matematic asociat mişcării studiate Sensul ei matematic este următorul: Cunoscându-se derivata doua funcţiei să se determine y Această cerinţă nu necesită în acest caz cunoştinţe speciale Luând succesiv de două ori primitiva ambilor membri ai ecuaţiei obţinem, rând pe rând , d d C C gt y t gt C t y = + = + Ultima expresie constituie soluţia generală ecuaţiei Observaţie Soluţia generală depinde în acest caz de două constante arbitrare, în timp ce în cazul ecuaţiei ea depindea doar de una O M y mg Ecuaţii diferenţiale ordinare de ordinul întâi IMPORTANT! Întotdeauna soluţia generală unei ecuaţii diferenţiale depinde de un număr de constante egal cu ordinul maxim de derivare al funcţiei necunoscute Vom reveni mai târziu cu justificări asupra acestui fapt semnificativ Să precizăm acum sensul fizic al constantelor şi Luând în prima expresie găsim d d v t y C t = = = unde este viteza iniţială particulei Analog, din doua expresie deducem = care reprezintă poziţia iniţială particulei Cu aceste noi notaţii pentru constante notaţii sugestive prin semnificaţia lor fizică soluţia generală ecuaţiei se pune sub forma v y gt y formă familiară cititorului încă din studiile liceale de fizică elementară Este clar acum care sunt datele suplimentare ce trebuie cunoscute pentru a determina acea soluţie care corespunde unei anumite mişcări, bine precizată: § poziţia iniţială particulei şi § viteza sa iniţială Se poate deci spune că satisface condiţiile d d , v t y y y = = Acestea se mai numesc şi condiţii iniţiale sau condiţii Cauchy Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Problema care constă în rezolvarea ecuaţiei astfel încât să satisfacă condiţiile iniţiale se numeşte problemă Cauchy sau problemă iniţială IMPORTANT! În cazul problemei Cauchy, condiţiile sunt puse în acelaşi punct! În exemplul de mai sus, în punctul Există însă situaţii în care acest tip de condiţii nu corespund fenomenului fizic Să luăm cazul unei bare simplu rezemate vezi figura de mai jos Problema constă în determinarea deflexiei încovoierii ca funcţie de x Nu vom intra în detalii de stabilire modelului matematic asociat Precizăm doar că acesta se prezintă sub forma ecuaţiei diferenţiale ordinare d d d d  + x y f x x y , numită şi ecuaţia Bernoulli-Euler O l y x Din figură se vede că la capetele şi ale barei deplasarea trebuie să fie nulă, adică y Condiţiile suplimentare se mai numesc şi condiţii bilocale Problema care constă în rezolvarea ecuaţiei cu condiţiile este o problemă bilocală sau problemă Picard engl : two-point problem Ecuaţii diferenţiale ordinare de ordinul întâi Aceste două tipuri de probleme asociate EDO sunt tipice şi acoperă mare parte din problemele mecanice şi fizice importante Din cele expuse mai sus se desprinde concluzia că nu se poate face un studiu sistematic de fenomen fizic fără se recurge la modelul său diferenţial După rezolvarea EDO sau EDP corespunzătoare, interpretarea soluţiei va permite cunoaşterea efectivă, previziunea şi deci controlul fenomenului studiat, iar acestea sunt deziderate majore ale ştiinţei FORMELE SUB CARE SE PREZINTĂ ECUAŢIILE DE ORDINUL ŞI SOLUŢIILE LOR Este evident faptul că ecuaţie diferenţială ordinară poate funcţiona doar în punctele în care este definită De exemplu, ecuaţia y are sens doar pentru Fiind dată ecuaţie diferenţială ordinară, trebuie determinat mai întâi domeniul pe care aceasta are sens; domeniul de definiţie al unei ecuaţii diferenţiale ordinare este cel al funcţiilor care definesc FORME ALE ECUAŢIILOR DE ORDINUL I A Forma generală ecuaţiilor diferenţiale ordinare de ordinul este, conform definiţiei şi relaţiei , F x y y d d unde este definit şi, de obicei, continuu în raport cu variabila independentă x, precum şi în raport cu funcţia necunoscută şi cu derivata acesteia, Forma generală se mai numeşte şi implicită, deoarece îl conţine implicit pe Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie B Dacă ¶ ¶ y F , atunci, conform teoremei funcţiilor implicite vezi cursul de Analiză Matematică, partea poate fi explicitat din şi obţinem forma canonică ecuaţiilor diferenţiale ordinare de ordinul I: y formă care se mai numeşte şi explicită C Dacă atunci se mai poate scrie y x , d d numită şi forma inversă, formă care poate fi folosită în vecinătatea acelor puncte ÎÂ în care tinde la infinit Evident, dacă nu tinde la infinit, formele şi sunt echivalente D Ecuaţia mai poate fi scrisă şi sub forma diferenţială: dy dx de asemenea echivalentă cu Forma diferenţială mai generală P dx dy este şi ea echivalentă cu fiecare dintre ecuaţiile Q P y x y , , d d P Q y y x , , d d ATENŢIE! În punctele , în care şi se anulează, nici una dintre ecuaţiile , nu este definită Ca şi în cazul ecuaţiei funcţiile şi sunt de cele mai multe ori continue pe domeniul de definiţie al ecuaţiei Ecuaţii diferenţiale ordinare de ordinul întâi E Forma simetrică EDO de ordinul este y X y x , d , d Fiecare din formele de mai sus pune în evidenţă anumite caracteristici şi posibilităţi de rezolvare ale ecuaţiilor de ordinul I Cel mai des întâlnite sunt formele şi FORME ALE SOLUŢIILOR Definiţia soluţie ecuaţiei diferenţiale în intervalul real [a,b] este o funcţie de clasă [a,b] care satisface identic adică y xÎ[a,b] Dacă există constantă astfel încât pentru orice xÎ[a,b], rezultă, evident, că este soluţie lui Ea se numeşte soluţie staţionară şi este deosebit de importantă pentru studiul calitativ al ecuaţiei Pentru rezolva ecuaţie diferenţială de ordinul se folosesc, după caz, formele menţionate în paragraful precedent şi, în funcţie de acestea, vom obţine şi soluţiile lor sub diferite forme Soluţiile unei EDO de ordinul pot fi determinate a sub formă explicită: xÎ[a,b]; b sub formă implicită: ; c sub formă parametrică: Î[ ]ÍÂ   = = t b y t x , , , Exemplu Funcţia y xÎ este soluţie explicită ecuaţiei Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie y x y VERIFICARE Într-adevăr, pe de parte d d d d x x x x x x x y - = - = = iar pe de altă parte, x y x - - Expresiile şi coincid Soluţia poate fi exprimată şi implicit: F VERIFICARE Într-adevăr, calculând diferenţiala lui găsim dF xdx ydy Din ultima egalitate deducem y x y adică tocmai Soluţia mai poate fi exprimată şi parametric: sin , cos , >   = = t y t x t VERIFICARE Putem scrie ecuaţia şi sub forma diferenţială xdx ydy Avem Ecuaţii diferenţiale ordinare de ordinul întâi d sin cos sin cos , d sin sin sin cos sin cos d d cos cos cos sin sin cos d x t y t x t + + + = = = de unde rezultă xdx ydy adică tocmai TIPURI DE ECUAŢII DIFERENŢIALE DE ORDINUL I REZOLVABILE PRIN CUADRATURI Există anumite ecuaţii de formă particulară, des întâlnite în aplicaţii, pentru care s-au găsit metode de rezolvare cu ajutorul cărora soluţia se exprimă folosind primitive ale unor funcţii Spunem, în acest caz, că ecuaţia se rezolvă prin cuadraturi integrări Vom aminti şi rezolva aici câteva asemenea tipuri de ecuaţii diferenţiale ordinare ECUAŢII CU VARIABILE SEPARATE Sunt de forma X dx dy unde şi sunt funcţii continue, depinzând de variabilele respectiv y MOD DE REZOLVARE Observăm că funcţia F ∫ dx ∫Y dy admite ca diferenţială membrul stâng al ecuaţiei Într-adevăr, y y F x x F dF d ¶ + ¶ ¶ = Rezultă deci dF astfel încât Prin urmare, soluţia generală a EDO este Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie ∫ dx ∫Y dy Exemplu Să se rezolve ecuaţia { { d e- = y x Y y X x x Rezolvare Este, evident ecuaţie cu variabile separate Calculând primitivele, găsim ln , d d , y y y Y y X x = = = ∫ ∫ ∫ deci soluţia generală EDO este - e- ln ECUAŢII CU VARIABILE SEPARABILE Acestea au forma P dx dy unde p,Q,q sunt funcţii continue în raport cu argumentele corespunzătoare MOD DE REZOLVARE Dacă pe domeniul de definiţie al ecuaţiei, împărţim cu şi obţinem dy q y Q y x p x P x care este ecuaţie cu variabile separate Conform cazului precedent, soluţia generală a ecuaţiei este Ecuaţii diferenţiale ordinare de ordinul întâi C q y Q y x p x P Exemplu Să se rezolve ecuaţia ydx dy Rezolvare Este EDO cu variabile separabile Împărţim cu şi, după simplificări, obţinem d d = - y y x x x Aceasta este ecuaţie cu variabile separate, deci soluţia generală este dată de y C y x x x = - ∫ d d sau, calculând primitivele, - ln valabilă pentru > I ECUAŢII DIFERENŢIALE OMOGENE, DE GRADUL m Definiţia funcţie :Â ®Â, se numeşte omogenă de gradul m dacă: f tx,ty tm "t ÎÂ Dacă egalitatea are loc doar pentru se numeşte pozitiv omogenă O ecuaţie omogenă de ordinul are forma P dx dy unde şi sunt omogene de acelaşi grad MOD DE REZOLVARE Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Facem schimbarea y xz®dy xdz zdx Introducând în ecuaţie, rezultă P xz dx xz xdz zdx Dar P,Qsunt omogene de gradul deci P xz xmP xz xmQ Rezultă xm[P dx xdz zdx ]= Mai departe, [P zQ ]dx xQ dz Aceasta este ecuaţie cu variabile separabile Împărţim cu zQ şi deducem , d = + + z P zQ z Q z x x deci soluţia generală ecuaţiei este, conform celor spuse mai sus, C P zQ z Q z x x = + ∫ ∫ d , d , Făcând notaţia ∫ + j z P zQ z Q z z d , , soluţia generală ecuaţiei se scrie astfel ln sau, trecând la exponenţială, Ecuaţii diferenţiale ordinare de ordinul întâi x Ce-j Revenind la variabila iniţială, obţinem soluţia generală EDO omogene sub forma j = x y x Ce Exemplu Să se determine soluţia generală pentru următoarea EDO omogenă: xy dx xy dy Rezolvare Avem xy iar xy Evident, această ecuaţie nu este nici cu variabile separate, nici separabile Să încercăm să verificăm dacă este omogenă, conform definiţiei : P xt, yt xy xy , Q tx,ty xy xy Deci ecuaţia este omogenă de gradul Pentru rezolvare, efectuăm schimbarea y xz®dy xdz zdx Rezultă succesiv [ ] [ ]d ; d , d , + = + = + = z z x x x x în final, obţinem zdx xdz care este ecuaţie cu variabile separabile Împărţind cu xz, găsim Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie d d = z z x x care este ecuaţie cu variabile separate Soluţia sa generală este z C x z x + d d , sau ln ln C Revenind la vechile variabile, avem ln ln y y x C x x + Trecând la exponenţială, rezultă soluţia generală ecuaţiei omogene C x y x x y × sau, altfel scris x y y Cx - ECUAŢII CU DIFERENŢIALE TOTALE EXACTE Sunt de forma P dx +Q dy Definiţia ecuaţie diferenţială ordinară de ordinul se numeşte ecuaţie cu diferenţiale totale exacte dacă există funcţie diferenţiabilă astfel încât dF dx dy Din Cursul de Analiză Matematică, partea I-a, se ştie că: Ecuaţii diferenţiale ordinare de ordinul întâi dF dx dy dacă şi numai dacă x Q y P ¶ = ¶ ¶ ¶ CONSECINŢĂ: Soluţia generală unei ecuaţii cu diferenţiale totale exacte este F unde este constantă arbitrară Deci rezolvarea unei ecuaţii cu diferenţiale totale exacte se reduce la determinarea unei funcţii de două variabile, atunci când se cunoaşte diferenţiala MOD DE REZOLVARE · Etapa Se calculează derivatele parţiale x Q y P ¶ ¶ ¶ ¶ ; dacă ele coincid, rezultă că ecuaţia este cu diferenţiale totale exacte, adică există astfel încât dF dx +Q dy · Etapa Deoarece diferenţiala unei funcţii este vezi Cursul de Analiză, partea y y F x x F dF d ¶ + ¶ ¶ = rezultă        = ¶ ¶ = ¶ ¶ , Q y F P x F Integrând prima relaţie în raport cu se obţine forma lui F: F x x = j Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie unde este funcţie arbitrară depinzând doar de y Derivând ambii membri ai acestei relaţii în raport cu vom avea y P y F x x + ¶ = ¶ ¶ ¶ , d unde este fixat, dar arbitrar ales, astfel încât , să aparţină domeniului pe care sunt definiţi şi Q Ţinând acum seama de condiţia deducem t Q y F x x + ¶ = ¶ ¶ ¶ , Comparând această relaţie cu expresia lui y F ¶ ¶ din rezultă Q de unde j şi deci expresia lui este j ∫ y y y t y fiind ales în aceleaşi condiţii ca În final, găsim pentru F F t y y x x , d = astfel încât soluţia generală ecuaţiei cu diferenţiale totale exacte se obţine sub forma Ecuaţii diferenţiale ordinare de ordinul întâi P C y y x x ∫ = unde este constantă arbitrară Dacă integrăm mai întâi doua relaţie în raport cu obţinem soluţia generală sub forma echivalentă cu P C y y x x ∫ , = Exemplu Să se determine soluţia generală pentru ecuaţia yexdx ex dy Rezolvare I Verificăm dacă este satisfăcută condiţia Avem    = + = x x Q y P y deci      = ¶ ¶ = ¶ ¶ e , , e , , x x x Q y y P y şi rezultă că ecuaţia este cu diferenţiale totale exacte II Aceasta înseamnă că există de clasă astfel încât        = + ¶ ¶ = ¶ ¶ e e , x x y y F y x F Din prima relaţie deducem Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie F yex de unde y F ¶ ¶ Egalând această expresie cu cea din rezultă ex ex de unde j y j Înlocuind această expresie în obţinem soluţia generală ecuaţiei : C y y = e Observaţii Acelaşi rezultat se obţine prin aplicarea directă formulelor generale de mai sus a Aplicăm formula în care se poate lua Obţinem , e e , d y y t y t y F t x t y t t x t t x y t x y + = + = = = = = = ∫ ∫ ∫ prin urmare soluţia generală ecuaţiei este tot C y y = e , Ecuaţii diferenţiale ordinare de ordinul întâi cu constantă arbitrară b Aplicăm acum formula şi aici se poate lua Obţinem y t x y x x t x y t F t = = = × ∫ + prin urmare soluţia generală ecuaţiei este aceeaşi C y y = e , cu constantă arbitrară FACTOR INTEGRANT Deoarece modul de rezolvare al unei ecuaţii cu diferenţiale totale exacte este extrem de simplu, s-au căutat căi pentru exploata şi în alte situaţii această idee extrem de atrăgătoare Fie ecuaţia P dx dy Ne putem pune următoarea PROBLEMĂ Dacă ecuaţia nu este cu diferenţiale totale exacte, am putea oare găsi funcţie , cu care, înmulţind-o, s-o transformăm într-o ecuaţie cu diferenţiale totale exacte? Funcţia se numeşte factor integrant Putem demonstra cu uşurinţă că: Există întotdeauna un factor integrant ecuaţie diferenţială ordinară de ordinul admite infinitate de factori integranţi Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Orice factor integrant al unei ecuaţii diferenţiale ordinare de ordinul este de forma unde este integrală sau, altfel spus, o soluţie ecuaţiei, iar este un factor integrant Dacă se cunosc doi factori integranţi ai unei ecuaţii diferenţiale ordinare de ordinul atunci soluţia acesteia se scrie fără cuadraturi CUM DETERMINĂM FACTORUL INTEGRANT? Presupunem problema rezolvată; am înmulţit deci ecuaţia cu funcţie m , obţinând mP dx mQ dy care este ecuaţie cu diferenţiale totale exacte Conform proprietăţilor diferenţialei vezi Cursul de Analiză, partea există , de clasă astfel încât dF mP dx mQ dy ceea ce implică Q y F P x F m ¶ = ¶ ¶ ¶ Dacă este de clasă atunci, evident, y x ¶ ¶ m ¶ ¶ deoarece derivatele sale mixte coincid, conform teoremei Schwartz vezi Cursul de Analiză, partea Derivând cele două produse, obţinem P Q P Q y x m ¶m ¶ ¶m ¶ ¶ ¶ Ecuaţii diferenţiale ordinare de ordinul întâi care este, de fapt, ecuaţie cu derivate parţiale, pe care trebuie s-o satisfacă factorul integrant m; am ajuns deci la problemă aparent mai complicată decât cea de la care am plecat Presupunem acum că unde este funcţie cunoscută ce depinde de şi Deoarece depinde de şi doar prin intermediul lui aplicăm regula derivării în lanţ: x ¶x × ¶w w = m ¶ ¶m d d , y ¶y × ¶w w = m ¶ ¶m d d Introducem aceste expresii în şi obţinem ¶ - ¶ ¶ = ¶ - ¶w ¶ ¶w w m y P x Q x Q y P d d sau ×m ¶ - ¶w ¶ ¶w ¶ - ¶ ¶ ¶ = w m j w x Q y P y P x Q d d Dacă noua expresie, notată este funcţie ce depinde doar de ecuaţia se scrie dm mdw Aceasta este ecuaţie cu variabile separabile, în şi w Împărţind cu deducem d d w m cu soluţia generală ln ∫j dw lnC deci Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie ∫ d C De fapt, ne interesează doar soluţie particulară ecuaţiei deci putem lua , de exemplu După ce am deteminat factorul integrant, înmulţim cu el ecuaţia dată şi obţinem o ecuaţie cu diferenţiale totale exacte, pe care rezolvăm conform modelului de la paragraful precedent Observaţie Acest mod de rezolvare depinde de alegerea funcţiei ω; alegerea depinde, la rândul ei, de abilitatea rezolvitorului Însă, de multe ori, are forme simple, sau este indicat Exemplu Să se rezolve ecuaţia d d = + + + + - y x y x y x y x P Q ştiind că admite un factor integrant de forma Rezolvare Calculăm y y P + = ¶ ¶ , Q x y ¶ ¶ + Rezultă că P Q y x ¶ ¶ ¶ ¶ astfel încât ecuaţia nu este cu diferenţiale totale exacte Căutăm un factor integrant de forma unde, conform indicaţiei, w Trebuie ca y x ¶ ¶ m ¶ ¶ Ecuaţii diferenţiale ordinare de ordinul întâi Avem d y d ¶m × ¶ w , d x d ¶m × ¶ w Pe de altă parte, din calculele de mai sus, rezultă Q P x y ¶ ¶ ¶ ¶ Din rezultă d Q P Q d x m ¶ × ¶ m w ¶ , deci x y d P P Q d y w - + - m ¶ ¶ Aceasta înseamnă că ×w= m w - m d d , care este ecuaţie cu variabile separabile Împărţind cu wm obţinem ecuaţia cu variabile separate w = w m dm d , pentru care, căutând soluţie particulară, găsim lnm lnw Rezultă că factorul integrant căutat este adică m Înmulţim deci ecuaţia cu Obţinem Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie p q  dx dy = Aceasta este ecuaţie cu diferenţiale totale exacte, căci p x y ¶ ¶ , q x x ¶ ¶ Căutăm funcţie astfel încât x y F q x y F p xy x F ⇒ j        = + ¶ ¶ = + ¶ ¶ Derivăm pe în raport cu : F xy x x ¶ j ¶ Trebuie deci ca - xy+j xy+ şi rezultă că j În final, funcţia are forma , y F ; soluţia generală ecuaţiei este deci y ECUAŢII DIFERENŢIALE LINEARE DE ORDINUL I Ecuaţia y p,qÎC ÍÂ Ecuaţii diferenţiale ordinare de ordinul întâi unde x y y d , defineşte ecuaţia diferenţială de ordinul lineară şi neomogenă Deci ecuaţia omogenă asociată este y A Membrul stâng al ecuaţiei defineşte operatorul care asociază fiecărei funcţii funcţia adică Ly De exemplu, dacă este definit ca Ly atunci el realizează următoarea corespondenţă de la funcţie la funcţie: e cos sin cos , , = ¾¾® = = ¾¾® + = ¾¾® = - x L L y Ly y Ly x y Ly x Putem spune că operatorul dat de este definit astfel: L :C ®C B Operatorul dat de este linear Definiţia Spunem că un operator :X®Y, unde sunt spaţii vectoriale reale/complexe, este linear dacă ax bx aL bL pentru orice ÎX şi orice reali/complecşi Dacă y , iar sunt constante reale/complexe, atunci Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie [ [ , def y y y y L y Ly Ly = b = + b a b deci este linear, conform definiţiei de mai sus Observaţie Recunoaştem un operator diferenţial linear după faptul că, întotdeauna în structura lui, atât funcţia necunoscută cât şi derivata ei sunt la puterea întâi Definiţia Numim nucleu al unui operator :X®Y şi notăm cu “ker” de la kernel, engl mulţimea elementelor din care îl anulează, adică ker {xÎX Y} Se ştie cursul de Algebră, anul că ker este subspaţiu vectorial al lui X Pentru operatorul diferenţial linear dat de evident ker {yÎC Ly } deci ker coincide cu mulţimea soluţiilor ecuaţiei lineare şi omogene Ecuaţia lineară şi omogenă poate fi scrisă sub forma unei ecuaţii cu variabile separabile: d d ⇒ = x y de unde, prin împărţire cu deducem succesiv ln d ln , d , d y c y x p x y y = + = = ∫ Ecuaţii diferenţiale ordinare de ordinul întâi În ultima expresie, este constantă arbitrară, pe care putem considera de forma lnC Trecând la exponenţială în ultima egalitate, rezultă y Ce-∫ dx care este soluţia generală ecuaţiei omogene asociate Observaţie Formula arată că dimensiunea subspaţiului vectorial ker L este În continuare, vom scrie ecuaţia sub forma Ly Putem demonstra imediat Teorema Soluţia generală ecuaţiei neomogene este suma dintre o soluţie particulară ecuaţiei neomogene şi soluţia generală ecuaţiei omogene asigurate Demonstraţie Într-adevăr, fie soluţie particulară ecuaţiei neomogene Aceasta înseamnă că LY Să efectuăm în schimbarea de funcţie y Introducând în obţinem Ly LY Lz Lz L = = + linear Dar Ly deci implică Lz adică zÎker ◘ CUM ÎL DETERMINĂM PE Y? Răspunsul la această întrebare îl dă Metoda variaţiei constantelor sau metoda lui Lagrange Căutăm pe de forma numai că va fi considerat funcţie de x: Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie ∫ x Y x d e Atunci ∫ ∫ x Y x d d e şi, înlocuind în ecuaţia neomogenă obţinem e e e d d d d ∫ ∫ ∫ ∫ - - + º p x p x p x LY x Însă LY deci C p ∫ = e ceea ce conduce la x C x d e deci se obţine prin integrare: C x p x e d d ∫ În final, soluţia particulară este obţinută direct prin cuadraturi Y x p x e d d ∫ ∫ Ţinând seama de teorema , rezultă că Soluţia generală ecuaţiei lineare şi neomogene se obţine prin cuadraturi şi este dată de y x p x e d d ∫ ∫ ∫ sau, echivalent, de Ecuaţii diferenţiale ordinare de ordinul întâi = ∫ ∫ ∫ y x p x e d d d unde este constantă arbitrară Pentru rezolvarea unei ecuaţii lineare de ordinul putem folosi deci una dintre ultimele două formule, însă în practică este mai simplu să procedăm direct Din cele spuse mai sus se desprinde următorul MOD DE REZOLVARE Etapa I Se asociază lui ecuaţia omogenă corespunzătoare: Lz Am arătat că soluţia generală acestei ecuaţii omogene este dată de formula deci z Ce-∫ dx Etapa II Conform teoremei , rămâne să determinăm pe soluţie particulară a ecuaţiei Aceasta se realizează cu metoda variaţiei constantelor, după cum am arătat Exemple Să se determine soluţia generală pentru următoarele ecuaţii: a xy Rezolvare Este ecuaţie diferenţială de ordinul lineară şi omogenă Ea se mai poate scrie succesiv d ; d d , , d d + = + = + = x x y y y xy x xy x y Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie ultima este ecuaţie cu variabile separate Soluţia ei generală este x C y y d ln d ∫ , sau C x y ln ln = , unde este constantă arbitrară Trecând la exponenţială, găsim e x y C - = b e x y xy x - Rezolvare Este ecuaţie diferenţială de ordinul lineară şi neomogenă Etapa Ecuaţia omogenă asociată este z xz Soluţia ei generală fost deja găsită la exemplul Ea este e x z C - = Etapa Pentru determina soluţie particulară ecuaţiei neomogene, folosim metoda variaţiei constantelor Căutăm pe de forma e x Y x - = Introducem în ecuaţia neomogenă: , e e e x x x x x Y xY xC xC x Y xC x x x - - - = + = şi cum Ecuaţii diferenţiale ordinare de ordinul întâi e x Y xY x - , rezultă că e e x x C x - , adică şi deci x C Obţinem e x x Y - = Soluţia generală ecuaţiei neomogene este aşadar e e x x x y C - = , sau, altfel scris e x x y C - = , unde este constantă arbitrară ECUAŢIA BERNOULLI Este de forma y ya, aÏ{ }, p,qÎC § Dacă rezultă ecuaţia de ordinul lineară şi neomogenă y § Dacă , rezultă ecuaţia de ordinul lineară şi omogenă py q Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie MOD DE REZOLVARE § Împărţim cu ya y p x y y × = a a- § Derivăm y y -a a- : a -a = a y y y y x d d Deci se transformă în y p x x y = - a- a- d d Notăm u ya- şi obţinem u a care este ecuaţie lineară şi neomogenă, având pe drept funcţie necunoscută O rezolvăm şi revenim la prin Exemplu Să se rezolve ecuaţia y y x y Recunoaştem în ea ecuaţie de tip Bernoulli, cu Rezolvare § Împărţim ecuaţia cu : Ecuaţii diferenţiale ordinare de ordinul întâi × = y y y Notăm y y u y u = Avem = x u care este ecuaţie lineară şi neomogenă § Rezolvăm ecuaţia lineară o Ecuaţia omogenă asociată este = x u Rezultă u x u = , de unde deducem ln lnC ln Soluţia generală ecuaţiei este u × o Soluţia generală ecuaţiei neomogene este deci u × +U unde este soluţie particulară lui pe care determinăm cu metoda variaţiei constantelor Rezultă, succesiv, x - U U x × = × - - = x U , Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie astfel încât C , adică x × - ; soluţia particulară este deci U Soluţia generală ecuaţiei este u × § Soluţia generală ecuaţiei Bernoulli este dată de y - = - × ECUAŢIA RICCATI Are forma y p,q,rÎC ÍÂ · Dacă rezultă ecuaţia lineară şi neomogenă y · Dacă rezultă ecuaţia Bernoulli Dacă se cunoaşte soluţie particulară , ecuaţia Riccati se rezolvă prin cuadraturi MOD DE REZOLVARE Într-adevăr, cu schimbarea Ecuaţii diferenţiale ordinare de ordinul întâi y avem şi, înlocuind în aceasta devine Y Yz Însă pY qY , deci satisface z [p Y]z care este ecuaţie Bernoulli, cu După rezolvarea ei, revenim la cu schimbarea de funcţie Exemplu Să se rezolve ecuaţia x y ştiind că admite soluţie particulară Y x x = Rezolvare Folosind schimbarea de funcţie y z x = obţinem x z x z x z x + Rezultă ecuaţia Bernoulli z z x z pe care am rezolvat-o la exemplul corespunzător cazului Bernoulli, găsind Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie z × Revenind la cu schimbarea rezultă soluţia generală ecuaţiei Riccati y x x - = - × , unde este constantă arbitrară Să menţionăm căteva cazuri particulare simple în care ecuaţia Riccati se rezolvă prin cuadraturi Dacă r xÎI atunci se arată că soluţia generală ecuaţiei Riccati este [ [ ∫[ ] = + j + j = ∫ ∫ x x C x C x y x d , e d d Presupunem, mai general, că b abq xÎI unde constantele şi nu sunt simultan nule Dacă atunci, cu schimbarea de funcţie y / obţinem pentru noua funcţie necunoscută ecuaţie Bernoulli u b a u  Ecuaţii diferenţiale ordinare de ordinul întâi Dacă şi sunt polinoame satisfăcând const atunci Y [P D] and [P D] sunt ambele soluţii ale ecuaţiei Riccati y Comentariu Ecuaţia Riccati este deosebit de importantă în aplicaţiile din mecanică, inginerie, fizică, chimie, etc ; de aceea, fost mult studiată Are serie de proprietăţi remarcabile de exemplu, oricare soluţii distincte ale unei ecuaţii Riccati date sunt totdeauna în raport anarmonic Sistemele de ecuaţii Riccati sunt printre cele mai des folosite în cercetări moderne din domeniul ştiinţelor naturii ECUAŢIA CLAIRAUT Este de forma y xy MOD DE REZOLVARE Folosim schimbarea p x y y = d d de unde rezultă imediat dy pdx Pe de altă parte, din rezultă y xp relaţie care, diferenţiată, devine Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie {y p p x d d d = Egalând cele două expresii ale lui dy, obţinem pdx= pdx+xdp dp deci dp Rezultă că cel puţin una din următoarele egalităţi este valabilă    + = = d , x p p o Cazul Dacă dp atunci şi deci y xC unde este constantă arbitrară Relaţia reprezintă soluţia generală a ecuaţiei Clairaut Geometric, soluţia ecuaţiei Clairaut reprezintă un fascicol de drepte o Cazul Dacă atunci şi deci y Rezultă    = j = , , y p x p care reprezintă ecuaţia parametrică unei curbe integrale pentru ecuaţia Clairaut, care nu se obţine din soluţia generală, particularizând pe C De aceea, această soluţie este o soluţie singulară Geometric, ea este înfăşurătoarea fascicolului de drepte reprezentat de soluţia generală Într-adevăr, dacă este un fascicol de curbe, atunci eliminând pe C între relaţiile Ecuaţii diferenţiale ordinare de ordinul întâi    = ¶ ¶ = x C C F F C obţinem înfăşurătoarea fascicolului În cazul ecuaţiei Clairaut, şi C F ¶ ¶ au următoarea formă    ¶ ¶ , x C C F F xC y Eliminând pe între cele două ecuaţii de mai sus, obţinem    = j = , , y C x C care sunt tocmai ecuaţiile parametrice ale soluţiei singulare Deci soluţia singulară ecuaţiei Clairaut este înfăşurătoarea fascicolului de drepte ce reprezintă soluţia sa generală Exemplu Să se rezolve ecuaţia y xy Rezolvare   = = ⇒ = = d d , p y x y p y p Egalând expresiile lui dy, deducem pdx+xdp pdp pdx⇒ dp , adică   = = d , x p p Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie o Cazul dp deci y ×C care reprezintă soluţia generală ecuaţiei Clairaut o Cazul Avem    = = = , , y x p de unde deducem imediat x y care reprezintă soluţia singulară ecuaţiei Clairaut O y x În figura de mai sus este înfăţişată soluţia singulară, tangentă în fiecare punct la una din dreptele fascicolului care reprezintă soluţia generală ecuaţiei Clairaut considerate ECUAŢIA LAGRANGE Este de forma A deci depinde linear de şi y Dacă împărţind cu el, obţinem Ecuaţii diferenţiale ordinare de ordinul întâi A C y g y A y B y y y = , Dacă , atunci este ecuaţie Clairaut; fost tratată în paragraful precedent Presupunem deci că MOD DE REZOLVARE Procedăm ca în cazul ecuaţiei Clairaut Fie deci p x y y = d d de unde rezultă imediat dy pdx Pe de altă parte, din rezultă y xf relaţie care, diferenţiată, devine { xf p p x d d d = Egalând cele două expresii ale lui dy, obţinem pdx= dx+xf dp dp deci [f p]dx+[xf ]dp Dacă const atunci ecuaţia este cu variabile separabile şi se rezolvă ca în paragraful În caz contrar, avem două situaţii posibile: a Atunci împărţim cu şi rezultă Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie d d = - + - + f p g p x f p f p p x Aceasta este ecuaţie diferenţială ordinară lineară şi neomogenă, cărei funcţie necunoscută este variabila independentă fiind p Rezolvând-o cu metoda descrisă la paragraful obţinem soluţia sub forma unde este o constantă arbitrară Din rezultă y [a ]f sau, altfel scris, y unde am folosit notaţiile a În final obţinem soluţia generală ecuaţiei Lagrange sub forma parametrică    = + = + , y p x p b Dacă admite soluţiile reale pi înlocuind în ecuaţia şi ţinând seama că pi pi rezultă soluţiile xpi pi relaţii care reprezintă ecuaţii ale unor drepte, pentru fiecare pi Aceste soluţii pot fi singulare Exemplu Să se determine soluţia generală ecuaţiei y Rezolvare Este ecuaţie de tip Lagrange Deci aplicăm schimbarea Ecuaţii diferenţiale ordinare de ordinul întâi y → dy pdx şi deducem y din care obţinem, prin diferenţiere, y dp d d Egalând cele două expresii ale lui dy, găsim p dp d d sau, după efectuarea calculelor, d  Rezultă că cel puţin una din următoarele egalităţi este valabilă:    = - = d , d p x p a Prima egalitate este de fapt ecuaţia cu variabile separate d dx , cu soluţia generală x C Din rezultă şi y C Soluţia generală ecuaţiei Lagrange se obţine deci în forma parametrică Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie        = + = + , y C x C Eliminând între cele două expresii din găsim soluţia generală sub forma implicită b Cea de doua egalitate implică , care, înlocuit în duce la soluţia singulară ecuaţiei Lagrange: METODA APROXIMAŢIILOR SUCCESIVE În paragraful precedent am pus în evidenţă unele tipuri de ecuaţii diferenţiale de ordinul care pot fi rezolvate prin cuadraturi, conducând la formule analitice concrete ale soluţiilor Dar nu sunt multe cazurile în care apar ecuaţii de aceste tipuri Acest neajuns ar putea fi compensat prin găsirea unor metode aproximative ale soluţiilor Una dintre cele mai uzitate asemenea metode este metoda aproximaţiilor succesive, sau metoda lui Picard Întrucât metoda este constructivă, vom prezenta în cadrul complet al teoremei de existenţă şi unicitate soluţiei problemei Cauchy TEOREMA CLASICĂ DE EXISTENŢĂ ŞI UNICITATE CAUCHYPICARD Teorema Fie problema Cauchy Ecuaţii diferenţiale ordinare de ordinul întâi    = = d d y f y x y Presupunem că satisface următoarele condiţii: ÎÂ b} C este Lipschitz în raport cu deci există constantă pozitivă astfel încât f £ " ÎW Atunci problema Cauchy admite soluţie unică yÎC unde este intervalul , lungimea sa fiind determinată astfel: M M b h a x y min sup , , ÎW =     = * Demonstraţie există şi este finit, căci este continuu pe compact Demonstrăm întâi EXISTENŢA SOLUŢIEI Integrând ecuaţia din şi ţinând cont de condiţia Cauchy, observăm că problema este echivalentă cu ecuaţia integrală ∫ x x y t Existenţa este constructivă, prin METODA APROXIMAŢIILOR SUCCESSIVE, care se mai numeşte şi metoda lui Picard, autorul ei Metoda este eficientă şi are un grad mare de aplicabilitate Ea poate fi utilizată şi în alte probleme Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie În cazul nostru, considerăm următorul şir aproximant pentru soluţia problemei : ÎN = + = + ∫ ∫ ∫ y n y t y t x x n n x x x x Urmăm câteva etape: I Demonstrăm că şirul } yn nÎN este bine definit şi toate funcţiile yn nÎN, au valorile numai în intervalul [y b], pentru orice xÎI II Arătăm că şirul } yn nÎN este uniform şi absolut convergent pe I În acest scop, considerăm seria S +K+ yn yn- +K ale cărei sume parţiale sunt chiar +K+ yn yn- yn Demonstrăm că §seria are termenii majoraţi de constante pozitive pe iar § seria numerică acestor constante este convergentă Conform criteriului lui Weierstrass vezi cursul de Analiză Matematică, partea I rezultă că Seria este absolut şi uniform convergentă pe I Să notăm suma acestei serii cu y Termenul general al lui este continuu, deci, conform proprietăţilor sumei seriilor de funcţii Cursul de Analiză Matematică, partea că Ecuaţii diferenţiale ordinare de ordinul întâi Suma seriei este continuă Se poate trece deci la limită în relaţia de definiţie şi avem ∫ x x y t ; cum şi sunt continue, rezultă că membrul drept al lui este derivabil, deci membrul stâng este de clasă În concluzie, satisface problema Cauchy UNICITATEA SOLUŢIEI Se demonstrează prin reducere la absurd ◘ PRINCIPIUL CONTRACŢIEI Metoda aproximaţiilor succesive aplicată ecuaţiilor diferenţiale ordinare implică un concept mult mai general, cu numeroase aplicaţii, anume, principiul contracţiei Îl vom prezenta pe scurt Fie mulţime pe care s-a definit distanţă metrică : d ´ ®Â cu proprietăţile: şi = , ,"x, yÎ proprietatea de simetrie, ,"x, zÎ inegalitatea triunghiului Astfel, împreună cu ce îndeplineşte proprietăţile de mai sus formează spaţiul metric Definiţii: Şirul {x nÎN este convergent în metrică către xÎ dacă şirul numeric { xn nÎN este convergent Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Şirul {x nÎN se numeşte Cauchy în metrică dacă pentru orice găsim un rang astfel încât d xn xn+ pentru orice rang şi orice pÎN se numeşte spaţiu metric complet dacă orice şir Cauchy în metrică admite limită în Să considerăm acum un operator ® Definiţii: se numeşte contracţie dacă există un număr pozitiv subunitar astfel încât d Tx,Ty rd "x, yÎ se numeşte punct fix pentru dacă x Tx Cu aceste precizări şi definiţii, putem enunţa acum, fără demonstra, Teorema Principiul contracţiei Fie un spaţiu metric complet şi T ® contracţie Atunci admite un punct fix unic APLICAŢIE: DEMONSTRAREA TEOREMEI CAUCHY-PICARD CU PRINCIPIUL CONTRACŢIEI Teorema Fie problema Cauchy    = = y y y Presupunem adevărate ipotezele teoremei , deci ÎC b} , ; Lipschitz în raport cu pe adică există constantă astfel încât Ecuaţii diferenţiale ordinare de ordinul întâi f atunci este funcţie descrescătoare şi lim t W t ®¥ a = deducem că volumul de lichid scade şi tinde spre valoare strict pozitivă această situaţie este prezentată în simularea de mai sus Dacă W a şi respectiv Figura Viteza rachetei în cazul frecării neliniare cu Figura Mişcarea rachetei în cazul frecării neliniare cu Figura Viteza rachetei în cazul frecării neliniare cu Figura Mişcarea rachetei în cazul frecării neliniare cu Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Figura Viteza rachetei în cazul frecării neliniare cu Figura Mişcarea rachetei în cazul frecării neliniare cu În toate simulările făcute am utilizat următoarele valori ale parametrilor şi condiţiei iniţiale: kg, kg ×m-a sa şi Aplicaţia M V Soare Problema fizică problema lui Cayley Să se studieze mişcarea unui unui corp solid de greutate care se deplasează pe un plan înclinat cu unghiul fiind legat de un lanţ înfăşurat fără frecare în figura Model matematic Aplicând teorema momentului, se obţine ecuaţia diferenţială de ordinul I X g p t v g P d d în care este masa totală sistemului mecanic la momentul fiind acceleraţia gravitaţiei, este acumularea de masă, este forţa exterioară, viteza la Ecuaţii diferenţiale ordinare de ordinul întâi momentul iar este viteza iniţială masei adiţionale Ecuaţia reprezintă modelul unui sistem mecanic de masă variabilă Fie greutatea lanţului pe unitatea de lungime; în acest caz, pentru orice deplasare greutăţii masa totală în mişcare va fi P qx Observăm că d d qv t P p Figura Sistem mecanic de masă variabilă Porţiunea înfăşurată lanţului fiind în repaus, putem considera că viteza iniţială a masei adiţionale este nulă Forţa exterioară este componenta după direcţia planului înclinat forţei aşa încât qx sin În felul acesta, ecuaţia capătă forma sin d d d d qx t P v t v P g Soluţie Ecuaţia care guvernează fenomenul este deci Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie sin sin d d Pv Pg qx a t Înmulţind la stânga cu Pv şi la dreapta cu qx dx dt şi apoi integrând, găsim sin qx C q g Pv Interpretare fizică Dacă admitem că pentru sistemul mecanic este în repaus la partea superioară planului înclinat, atunci, din condiţia rezultă că = sin a C iar pătratul vitezei este dat de sin sin a + a + + = P qx gx qx x P qx P qx P q g v În cazul particular lanţul cade liber se obţine sin , d d gx t x v de unde sin ; d t g x x a apoi sinat , aşa încât sin , sin , sin , = g t t g t t g x mişcarea elementelor lanţului fiind uniform accelerată Ecuaţii diferenţiale ordinare de ordinul întâi Aplicaţia M V Soare Problema fizică Să se determine forma de echilibru unui fir elastic suspendat între două puncte, având aria secţiunii şi modulul de elasticitate al materialului E Va fi studiat cazul greutăţii proprii firului mg ] Model matematic Fie efortul în fir şi componentele sale după axele Ox şi Oy: S dx ds dy ds respectiv figura În starea deformată firului, ecuaţiile de proiecţie pe cele două axe dau d d d d s x S s , d d d d g EA S s y S s = unde este greutatea proprie firului pe unitatea de lungime se ia masa egală cu unitatea Din rezultă x s S S s x S d d const, d d = şi, introducând în obţinem d d d d d d g x s EA S s y s S + Ţinând seama că ds dx dy dx rezultă ecuaţia diferenţială nelineară de ordinul I d d g EA y S x y S + + Soluţie Notăm şi considerăm ca variabilă independentă; obţinem d d + = + EA p S g S p x Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Figura Deformarea unui fir elastic suspendat între două puncte Prin integrare rezultă ln  + C EA S g S x Deoarece pentru avem , deducem şi ln  + EA S g S x Dacă înmulţim cu dy dx rezultă d d + = + p p p EA S g S p y Prin integrare, deducem C EA S g S y Deoarece pentru avem , rezultă şi deci EA S g S y Relaţiile şi constituie representarea parametrică fibrei deformate Se vede că atunci când EA®¥ cazul firului inextensibil se regăseşte ecuaţia lănţişorului Ecuaţii diferenţiale ordinare de ordinul întâi Efortul poate fi determinat dintr-o condiţie geometrică legată de lungimea totală firului Aplicaţia M V Soare Problemă Considerăm starea de eforturi de membrană simetrică intr-o placă subţire de rotaţie, supusă la încărcare exterioară de componente după tangenta la meridian, respectiv după normala la suprafaţa mediană şe cer expresiile generale ale eforturilor meridiane şi inelare Nj Nq figura Model matematic Ecuaţiile de echilibru unui element de placă sunt cos , d d = j Yr j Z r N r N Variabila independentă problemei este unghiul meridian măsurat în sens direct orar de la creştet, fiind unghiul de-a lungul cercului paralel Alte mărimi fizice implicate în model sunt: raza cercului paralel, raza de curbură r cosj dr dj curbei meridiane prima rază de curbură principală suprafeţei mediane şi sin doua rază de curbură principală suprafeţei mediane Soluţie Deoarece ecuaţia este algebrică, avem de fapt singură ecuaţie diferenţială în Nj pe care obţinem eliminând pe Nq determinat din cu expresia Zr r r Nq Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Figura Eforturile de membrană într-o placă subţire de rotaţie Introducând în deducem cos d d = j Yr Zr r Ţinând seama de relaţiile dintre razele şi rezultă sin sin cos d d N j j Notând acum Njr sin obţinem pentru ecuaţie diferenţială ordinară de ordinul lineară şi neomogenă, studiată la sin cos d d Y r y j j Integrând-o, rezultă C r N + j j sin cos d sin , C fiind constantă arbitrară Efortul inelar se obţine direct din sin cos d sin r r Y C r N Zr j q ∫ Ecuaţii diferenţiale ordinare de ordinul întâi Constanta poate fi determinată dintr-o condiţie impusă la marginea superioară js sau la creştet Aplicaţia M V Soare Problemă Se cere să se determine starea de tensiune normală, ca funcţie de timp, pentru un corp Maxwell Model matematic Pentru explicarea relaxării se alcătuieşte modelul Maxwell, prin combinarea unui model Hooke elastic şi unui model Newton vâscos figura Starea de tensiune rezultă ca sumă stărilor de deformaţie celor două corpuri; astfel, tensiunea totală const este compusă din § deformaţia elastică arcului, dată de eelastic / unde este modulul longitudinal de elasticitate, şi din § deformaţia vâscoasă, evascos Prin urmare figura evascos s E Derivând în raport cu timpul e& , obţinem s evascos & & E Ştiind că pentru corpul Newton subzistă relaţia vascos h e & în care prin s-a notat coeficientul de vâscozitate dinamică, care este constant Astfel, devine s h s E & Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Soluţia şi interpretarea ei fizică Ecuaţia diferenţială este lineară şi omogenă, adică de tipul celor studiate în Separând variabilele, obţinem t E d d h = s s , ceea ce implică ln ln t, E C h s unde este constantă arbitrară Soluţia generală ecuaţiei este deci e t E C h - s = Presupunem îndeplinită condiţia iniţială s Rezultă Soluţia problemei Cauchy considerate este t R h - s Figura Modelul Maxwell Variaţia lui în funcţie de t Variaţia lui ca funcţie de este dată în figura Graficul reprezintă o exponenţială descrescătoare care admite ca asimptotă axa timpului Ecuaţii diferenţiale ordinare de ordinul întâi Aplicaţia M V Soare Problema fizică Un fir trece peste un scripete circular fix, de rază figura între fir şi scripete luând naştere forţă de coeficient de frecare de alunecare f Dacă la una din extremităţile firului acţionează tensiune ce tensiune trebuie să se exercite la cealaltă extremitate pentru ca firul să înceapă să alunece pe scripete? Model matematic deoarece scripetele este rugos, reacţiunea ds asupra unui element de fir va avea pe lângă componentă normală ds şi una tangenţială Φ ds numită forţă de frecare de alunecare Din echilibrul unui element de fir figura se obţine ecuaţia vectorială dT+ ds mai putem scrie Tτ Nν fNτ ds d Figura Echilibrul unui fir pe un scripete În este reacţiunea normală de-a lungul vectorului unitate iar fN este reacţiunea tangenţială la limită de-a lungul vectorului unitate τ Figura Eforturile acţionând pe arcul s Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie În definitiv, utilizând formula lui Frenet dτ ds / putem scrie sistemul de ecuaţii care modelează fenomenul , d d - = - = N R T fN s T Soluţie Eliminând reacţiunea normală deducem următoarea ecuaţie diferenţială prdinară de ordinul lineară şi omogenă ds Rdq d d = q q fT T Conform ipotezei, ecuaţia trebuie integrată cu condiţia iniţială Cu metoda de la obţinem imediat soluţia generală ecuaţiei sub forma T Ce fq în care este constantă arbitrară Condiţia iniţială conduce la determinarea soluţiei problemei Cauchy , T fq Interpretare fizică Pentru putem scrie fa unde s-au pus în evidenţă mărimile tensiunilor la capetele firului Echilibrul poate avea loc şi pentru T ; în acest caz, forţa de frecare de alunecare îşi schimbă sensul şi rezultă T fa Se obţine astfel condiţia de echilibru lui Euler - fa Î = for , for t T Q T Qe unde şi sunt constante Model matematic Pentru deduce ecuaţia diferenţială care guvernează mişcarea, observăm că, în intervalul de timp dt, suma dintre volumul acumulat şi cel evacuat este egală cu volumul afluent Adh Ch dt Qedt Figura Chiuveta unui lac de acumulare Soluţie Pentru primul interval vom scrie ecuaţia sub forma t Q Ch A h e d d = - Introducând notaţia Qe schimbarea de funcţie h dh ydy conduce la ecuaţia cu variabile separate vezi Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie t y y y C A d d = b Descompunând fracţia precedentă în fracţii simple , b + - b b + + b b = b + + b = b y y y y y y y y y y ecuaţia diferenţială devine y t y y y C y A d d b + - b b + + b + b Integrând, obţinem ln arctan ln t t y y y C A  b - + b b , unde este constantă de integrare Soluţia precedentă se mai scrie ln arctan t t y y y y C A b - b b + b b ; revenind la funcţia iniţială găsim ln arctan t t h h h h C A b - b b + b b Pentru avem deci arctan t C A b - În final, obţinem Ecuaţii diferenţiale ordinare de ordinul întâi t T] h h h h h C A ln arctan = + b - b + b b La momentul devine ecuaţie transcendentă T h h h h h C A T T T T = + b - b + b b ln arctan care determină nivelul hT al apei Pentru avem Qe şi ecuaţia ia forma mai simplă Adh +Ch dt sau h- dh dt C A Integrând, obţinem h t C A în care este constantă de integrare; ea se determină din condiţia hT Rezultă h t C A T Astfel, deducem soluţia finală T C A t , , unde t T t T A C h h T  = - , Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Aplicaţia M V Soare Problema fizică Un recipient având aria secţiunii transversale orizontale are pe fund un orificiu care poate evacua un debit Qd Ch unde este constantă iar h este adâncimea apei din recipient Se cere să se studieze variaţia în timp nivelului al apei din recipient dacă debitul afluent Qe se prezintă în următoarele variante iniţial vasul este gol, adică pentru avem : a [ ]   > Î = for , for e T Q T Q b          Î = , , for , for , e T t T t Q T t T t Q Q unde şi sunt constante În figura se dă schema de calcul, iar în figura sunt date cele două legi de variaţie lui Qe Model matematic Pentru obţine ecuaţia diferenţială care guvernează mişcarea, observăm că, într-un interval de timp dt, suma dintre volumul acumulat şi cel evacuat este egală cu volumul afluent e d d Ch Q t h A Aceasta este ecuaţie diferenţială de ordinul nelineară şi neomogenă Soluţie Cu schimbarea de funcţie h dh ydy ecuaţia devine Ecuaţii diferenţiale ordinare de ordinul întâi e d d Cy Q t y Ay şi putem trece la examinarea cazurilor din enunţ a Pentru Î[ ,T ecuaţia este cu variabile separate vezi Q Cy Ay y t - = d d Introducând notaţia soluţia generală ecuaţiei precedente devine + ln t A C y , unde este constantă de integrare; revenind la variabila soluţia precedentă se scrie ln + t A C h C Q C Q h Introducând condiţia iniţială pentru rezultă C Q C Q C A ln , astfel încât devine t ] A C Q Ch C Q h , ln + În particular, la momentul , avem T A C Q Ch C Q h T T ln + ; relaţia determină înălţimea hT Pentru intervalul , ecuaţia diferenţială ia forma Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie d d Ch t h A sau d d = h h A , cu soluţia generală Ah Ct în care este constantă de integrare, ce urmează fi determinată din condiţia de continuitate; pentru trebuie să avem hT determinat de relaţia Figura Recipient cu orificiu: schema de calcul; legile de variaţie ale lui Qe Ah CT Introducând expresia lui în obţinem Ah Ct Ah CT , care determină explicit nivelul h T A C h hT  , Ecuaţii diferenţiale ordinare de ordinul întâi Timpul poate fi explicitat din şi sub forma [ ]        + > Î  - = pentru ln pentru h T C A T t T Q Ch C Q h C A t T b Ecuaţia diferenţială devine T t Ch Q t h A d d + , pentru primul interval; cu schimbarea de funcţie ecuaţia se mai scrie T t A tu Ct Q + Simplificând cu rezultă ecuaţia cu variaile separate t t C Au T Q A d d = - , deci K K C Au T Q A u t ln ln d ln + + - = , unde este constantă pozitivă arbitrară Primitiva din membrul drept se scrie astfel ln ln , d d d d v v v v v v v v v v v v v v v v v v A v Av v T Q Av Cv Av v F u - - - - = - - - = - = + = ∫ ∫ ∫ ∫ Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie în care şi ,v sunt rădăcinile ecuaţiei algebrice = T Q Av Cv , şi sunt întotdeauna reale În consecinţă, , obţinem d d d d d d e d d , d d d d e d d , d d e d d = + = = - - - t y t y t y x y t y t y x y t y x y t t t deci Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie d d d d d d d d , d d d d d d , d d d d t y t y t y x y x t y t y x y x t y x y x = + = = Prin urmare k k k x y x d d se exprimă numai cu t y d d , , k k t y d d , iar ecuaţia se transformă în , d d d d , d d , = - t y t y t y F Luând p t y = d d şi ca variabilă independentă, obţinem ecuaţie având ordinul redus cu unitate Exemplu Pentru studia înconvoierea unei plăci subţiri, circulare de rază încastrate pe contur şi supusă unei sarcini concentrate în centrul ei, se utilizează ecuaţia diferenţială y kx R] x y x x y x d d d , unde este constantă Modelul matematic este reprezentat de ecuaţie diferenţială ordinară lineară şi neomogenă, de ordinul II Este chiar ecuaţie de tip Euler Pentru rezolva ecuaţia omogenă asociată facem schimbarea de variabilă et procedând la fel ca mai sus; de altfel, acesta este şi modul de rezolvare al ecuaţiei Euler vezi paragraful Obţinem d d - = t y , Ecuaţii diferenţiale ordinare de ordinul n care are soluţia , x x C y unde sunt constante reale Altă metodă de rezolvare ecuaţiei omogene este utilizarea substituţiei x y u , care permite reducerea ordinului cu unitate Obţinem xu''+ u'= Notând u'= p, obţinem ecuaţia diferenţială de ordinal întâi cu variabile separabile xp'+ Rezultă , x x C C y x C u x D p = + unde , K C şi sunt constante reale Soluţia generală ecuaţiei neomogene se află utilizând metoda variaţiei constantelor şi este x x k x x C k y ln ECUAŢII REDUCTIBILE LA EDO CU COEFICIENŢI CONSTANŢI Dacă, prin intermediul unei schimbări de variabilă sau funcţie, reuşim să transformăm ecuaţie diferenţială ordinară într-una lineară şi cu coeficienţi constanţi, atunci, prin transformarea inversă, putem exprima soluţia ecuaţiei date pornind de la cea a ecuaţiei transformate, pe care ştim să rezolvăm Vom da câteva exemple edificatoare Fie ecuaţia d d d d x y x x y Ly Să efectuăm schimbarea de variabilă Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie x cost Reconstituim ecuaţia: , d d sin cos d d sin cos d d sin cos d d sin d d sin cos d d sin d d sin d d d d sin d d d d d d = = ´ ´ ´ + = = = ´ ´ ´ n y t y t t t y t t t y Ly t t n t y t t y t t t y x t y t y x t t t y x y y y x x n şi rezultă d d + = t y Aceasta este ecuaţie diferenţial ordinară lineară şi omogenă, cu coeficienţi constanţi Conform celor arătate anterior vezi şi aplicaţia despre oscilatorul armonic de la paragraful un sistem fundamental de soluţii este y cosnt, sin nt sau, revenind la variabila x, y cosn arccos sin arccos Dacă , atunci cos arccos Pornind de la această observaţie, se demonstrează că pentru impar, este polinom de gradul în x Aceste polinoame sunt polinoamele Cebîşev Fie ecuaţia Bessel, mult folosită în aplicaţii inginereşti: x xy Ecuaţii diferenţiale ordinare de ordinul n Este ecuaţie lineară, de ordinul II, cu coeficienţi variabili Soluţia sa se caută sub formă de serie şi are drept rezultat introducerea funcţiilor Bessel, care depind de indicele ν Să scriem ecuaţia Bessel pentru n : Ly xy Aplicăm acestei ecuaţii schimbarea de funcţie x z y Reconstituim ecuaţia : , = + = ´ ´ × + = - - - - Ly z x x x y z y z y z de unde deducem pentru ecuaţia diferenţială lineară şi omogenă, cu coeficienţi constanţi, z Soluţia generală acestei ecuaţii este z cos sin revenind la prin obţinem Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie x c x y = cos sin , care este soluţia generală ecuaţiei Bessel pentru indicele ECUAŢIA EULER Fie ecuaţia diferenţială ordinară lineară, de ordinul cu coeficienţi variabili: = a xy n n Aceasta este ecuaţia Euler; observăm că derivatele de ordinul ale lui sunt înmulţite cu puteri de acelaşi ordin ale lui x Îi vom aplica schimbarea de variabilă x et Pentru mai bună înţelegere, vom face acest calcul pentru cazul cazul n arbitrar se tratează absolut analog Fie deci EDO a xy Avem t y x t t y x y t d d e d d d d d d × → x t t d d e d d Reconstituim ecuaţia, notând cu operatorul identitate Ey : d d d d d d d d d d d d e d d d d d d e d d d d e d d e d d e d d d d e d d e d d e d d e d d e d d =  = ´ = ´ = - - - y y t E a t t E a t E t t Ly a E a t E t t y t y x t y a x E a t t y x t y a x a t y x y a x a a t t t t t t Ecuaţii diferenţiale ordinare de ordinul n Rezultă EDO lineară, cu coeficienţi constanţi d d d d d d d d d d d d y t E a t t E a t E t t Ly După cum am arătat în paragraful , căutăm soluţii de forma y ert Ecuaţia caracteristică asociată este a După ce rezolvăm, găsim un sistem fundamental de soluţii şi scriem soluţia generală EDO Observaţie Combinând schimbarea de variabilă cu forma exponenţială , constatăm că y rt xr r ln eln Deci în aplicaţii este mai simplu să căutăm direct soluţii de forma y xr Exemplu Să se determine soluţia generală ecuaţiei x xy Soluţie Este ecuaţie de tip Euler şi deci căutăm direct soluţii de forma Avem [ ] , = = ´ ´ = ´ = - - r r r r Ly x x x x rx y x deci ecuaţia caracteristică este Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie r şi are rădăcinile Soluţiile corespunzătoare vor fi i i x- x- Pentru rămâne în cadrul real, folosim formulele lui Euler Putem scrie x y x x y x x ln x i ln i e , e = = astfel încât, luând partea reală şi imaginară lui obţinem sistemul fundamental real x x Y x x Y sin ln , cos ln ; soluţia generală ecuaţiei Euler este x c x y x = cos ln sin ln APLICAŢII ÎN MECANICĂ, FIZICĂ ŞI INGINERIE Aplicaţia Oscilaţii liniare D Comănescu, I Caşu Problema fizică Considerăm un corp material supus unei forţe de tip elastic Fe ® Vom privi corpul ca un punct material cu masa constantă m Pe baza cunoscutei legi a lui Hooke forţa elastică este direct proporţională cu vectorul de mişcare x r Modelul matematic al mişcărilor provine din legea II-a lui Newton În aplicaţiile practice intervin şi alte forţe asupra punctului material Atunci când acesta se mişcă printr-un mediu rezistent semnificativă este şi forţa de frecare, pe care noi vom considera direct proporţională cu viteza punctului material, mai exact F v ® ® = × Constanta de frecare este strict pozitivă din considerente fizice Dacă Ecuaţii diferenţiale ordinare de ordinul n punctul material se află în vecinătatea suprafeţei Pământului atunci trebuie să ţinem seama şi de greutatea punctului material situaţie foarte generală este aceea în care punctul material se află şi într-un câmp exterior variabil ® În cele ce urmează studiem mai multe situaţii în care forţa rezultantă cuprinde întotdeauna forţa elastică În această aplicaţie ne vom ocupa cu studiul mişcărilor rectilinii ale punctului material Notăm cu componenta mişcării pe dreapta aleasă Oscilaţii libere În această secţiune considerăm că forţa elastică este singura forţă ce acţionează asupra punctului material Modelul matematic al mişcărilor liniare este: m x ·· × , unde este constantă strict pozitivă numită constanta elastică Ecuaţia diferenţială poate fi scrisă sub forma m ·· × × care este ecuaţie diferenţială de ordinul liniară şi omogenă vezi Ecuaţia caracteristică asociată este m×l Rădăcinile acesteia sunt k i m l În fizică expresia k m se notează cu şi este numită frecvenţa unghiulară oscilaţiei sau pe scurt frecvenţă Soluţia generală a ecuaţiei diferenţiale este cos ×t sin ×t Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie unde şi sunt constante reale arbitrare Aceste constante pot fi unic determinate cunoscând poziţia iniţială şi viteza iniţială Mişcarea determinată de aceste condiţii iniţiale este soluţia unică următoarei probleme Cauchy: x x x x x v w ·· ·  ×     = Această soluţie are expresia cos sin v x t w = × Figura Oscilaţii libere Analizând expresia mişcării punctului material se observă urmatoarele proprietăţi:  mişcarea este mărginită iar maximul funcţiei de mişcare se numeşte amplitudinea oscilaţiei şi are valoarea v x w + ; Ecuaţii diferenţiale ordinare de ordinul n  mişcarea este periodică de perioadă T p w = În Figura este prezentată simulare numerică, făcută cu programul MAPLE mişcării punctului material Oscilaţii amortizate Punctul material este acţionat de forţa elastică şi forţă de frecare cu un mediu rezistent Modelul matematic al mişcărilor liniare este: m x ·· · × Ecuaţia diferenţială poate fi scrisă sub forma m ·· · × × × care este ecuaţie diferenţială de ordinul liniară şi omogenă Ecuaţia caracteristică asociată este m×l +m ×l şi are rădăcinile , m k m m m m l l - × × × = = × Deosebim trei cazuri: ×m× În acest caz l şi l sunt valori reale strict negative Soluţia generală ecuaţiei diferenţiale este exp ×t exp ×t Observăm că toate soluţiile tind spre când ®¥ Figura prezintă o simulare mişcării ii ×m× În acest caz soluţia generală ecuaţiei diferenţiale este Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie exp t x t m = × × × Şi în acest caz toate soluţiile tind spre când ®¥ Figura Oscilaţii amortizate Cazul ×m× Figura Oscilaţii amortizate Cazul ×m× Ecuaţii diferenţiale ordinare de ordinul n iii × În acest caz l şi l sunt complex conjugate Soluţia generală ecuaţiei diferenţiale este exp cos sin t x t m = × × × × unde am notat prin m k m m w - × × = × Toate soluţiile tind spre când ®¥ Figura Oscilaţii amortizate Cazul ×m× Oscilaţii forţate neamortizate În situaţia de faţă, asupra punctului material acţionează forţa elastică şi un câmp exterior variabil ® Presupunem că acest câmp este dirijat de-a lungul dreptei pe care se mişcă punctul material şi notăm cu componenta corespunzătoare Modelul matematic al mişcărilor liniare este: m ·· × Ecuaţia diferenţială poate fi scrisă sub forma Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie m ·· × × care este ecuaţie diferenţială de ordinul liniară şi neomogenă Vom studia această ecuaţie în câteva situaţii importante din punct de vedere practic Vom presupune că la momentul iniţial punctul material este în repaus şi în echilibru ceea ce din punct de vedere matematic înseamnă , x · = m× constant Această situaţie corespunde cazului în care punctul material se află sub influenţa forţei elastice şi greutăţii provocate de Pământ Constanta este acceleraţia gravitaţională Problema Cauchy mişcării este x mg x x w ·· ·  ×     = unde am utilizat notaţia făcută în paragraful referitor la oscilaţii libere Soluţia acestei probleme este cos g x t w = Remarcăm că punctul material execută mişcare mărginită şi periodică, de perioadă T p w = ×t cu constantă Această situaţie poate fi întâlnită atunci când aproximăm câmpul exterior de forţă prin partea sa liniară Problema Cauchy mişcării este Ecuaţii diferenţiale ordinare de ordinul n a x m x x w ·· ·  × ×    =   = Soluţia acestei probleme este sin a x t m w w w = × × × Se observă că oscilaţia este nemărginită × exp ×t cu şi constante Mişcării îi corespunde problema Cauchy exp a x m x x w a ·· ·  × × ×    =   = Este ecuaţie diferenţială ordinară de ordinul II, cu coeficienţi constanţi, lineară şi omogenă Soluţia acestei probleme este exp cos sin a x t m a w w w = × × Cazul de rezonanţă × cos ×t cu constantă În acest caz frecvenţa câmpului exterior coincide cu frecvenţa unghiulară a oscilaţiei se vedea secţiunea de oscilaţii libere Mişcarea este modelată de problema Cauchy cos a x m x x w w ·· ·  × ×    =   = Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Ecuaţia diferenţială este de ordinul II, lineară şi neomogenă, cu coeficienţi constanţi Soluţia acestei probleme este sin a x t m w w = × × × În figura este prezentată simulare numerică mişcării punctului material din care se observă că amplitudinea oscilaţiei creşte, mişcarea fiind nemărginită Acest fenomen este responsabil pentru numeroase catastrofe tehnice Figura Oscilaţii forţate, neamortizate cu rezonanţă Cazul de non-rezonanţă × cos ×t cu constante şi w În acest caz frecvenţa câmpului exterior nu coincide cu frecvenţa unghiulară oscilaţiei Modelul mişcării este reprezentat de problema Cauchy cos a x m x x w s ·· ·  × ×    =   = asemănător cazurilor precedente Ecuaţii diferenţiale ordinare de ordinul n Figura Oscilaţii forţate, neamortizate, cazul de non-rezonanţă Soluţia acestei probleme este cos cos a x t m w s w s = × × Observăm că mişcarea punctului material este mărginită şi este suma două funcţii periodice de perioade diferite În figura este prezentată simulare numerică a mişcării punctului material Aplicaţia Mişcarea pendulului simplu G Cosovici, S Comşa Pendulul simplu este un punct material cu masa suspendat de articulaţie fixă O prin intermediul unui fir inextensibil şi fără greutate având lungimea figura La momentul pendulul se află într-o configuraţie de repaus, în care firul formează cu verticala unghiul > Din această poziţie, el este lăsat să se mişte liber În afară de greutatea proprie, asupra punctului material acţionează tensiunea din fir vom neglija frânarea exercitată de aerul atmosferic Aceste forţe definesc planul traiectoriei parcurse de pendul Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Figura Modelul mecanic al unui pendul simplu Model matematic Dacă se utilizează un sistem de coordonate polare centrat în articulaţia figura mişcarea punctului material poate fi descrisă cu ajutorul funcţiei q =q care defineşte unghiul curent format de fir cu verticala pe tot parcursul discuţiei care urmează, vom considera că unghiul este măsurat în radiani În raport cu acest reper, singurele deplasări ale pendulului au loc pe direcţie circumferenţială, direcţie în lungul căreia acţionează componentă greutăţii neechilibrată de tensiunea din fir Observaţia de mai sus ne permite să scriem următoarea expresie celui de al doilea principiu al dinamicii figura : mlq&&+ mg sinq În egalitatea este acceleraţia gravitaţională După simplificarea cu m, devine O x q mg mg cosq mg sinq T mg cosq q q l Ecuaţii diferenţiale ordinare de ordinul n lq&&+ sinq Datorită neliniarităţii, ecuaţia diferenţială nu poate fi rezolvată analitic Totuşi, pentru unghiuri mici maximum rad , dezvoltarea în serie sin ! q =q +q ne permite să operăm cu aproximarea sinq »q Cu ajutorul acesteia, ecuaţia diferenţială se rescrie sub forma liniarizată g q&&+ = l Soluţie Ecuaţia este ecuaţie diferenţială lineară şi omogenă, de ordinul II, cu coeficienţi constanţi Soluţia sa generală are expresia vezi cos , g q t = în care > şi j sunt constante Valorile lor se determină impunând condiţiile iniţiale =q q& Particularizarea funcţiei =q definită sub forma pentru cele două constrângeri de mai sus, conduce la sistemul de ecuaţii Acosj =q lAsinj Prin rezolvarea acestuia se obţine =q Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie După înlocuirea expresiilor ale constantelor şi formula devine =q cos Interpretare fizică Se observă că soluţia descrie evoluţie periodică a unghiului =q Valorile extreme pe care le ia funcţia =q sunt Timpul care separă două treceri succesive printr-un maximum sau minimum se numeşte perioadă şi este calculabil cu formula T După cum se poate remarca, este constantă Oscilaţiile pendulului sunt aşadar izocrone Această concluzie este totuşi valabilă numai în ipoteza oscilaţiilor de amplitudine foarte mică Aplicaţia Mişcarea pendulului simplu în prezenţa frânării exercitate de aerul atmosferic G Cosovici, S Comşa Vom relua exemplul precedent, ţinând cont şi de rezistenţa aerului atmosferic Admitem că frânarea este proporţională cu viteza punctului material De asemenea, adoptăm ipoteza micilor oscilaţii Model matematic În atare condiţii, principiul al doilea al dinamicii se scrie sub forma mlq&&+ bq& mgq unde este rezistenţa aerului presupusă constantă Pentru comoditatea calculelor, este convenabilă definirea coeficientului de frânare b m g × l Cu ajutorul acestei mărimi, ecuaţia diferenţială devine Ecuaţii diferenţiale ordinare de ordinul n q&&+ q& = l Soluţie Ecuaţia de ordinul II este lineară şi omogenă, cu coeficienţi constanţi Soluţia sa generală are expresia vezi Ael el unde şi rezultă din condiţiile iar l şi l sunt soluţii ale ecuaţiei caracteristice l + deci Pendulul va efectua mişcare periodică numai atunci când g Pentru acest caz, soluţiile se pot rescrie sub forma iw, iw unde g g l Soluţia generală devine atunci e-g Ae-iwt eiwt Întrucât putem găsi oricând două constante şi j care să garanteze satisfacerea egalităţilor A i A e- ej Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie relaţia admite rescrierea sub forma echivalentă Ae-g cos wt +j Aplicând constrângerile soluţiei generale deducem sistemul de ecuaţii Acosj =q Acosj Asinj Prin rezolvarea acestuia obţinem cos arctg arctg q g w w × După înlocuirea expresiilor ale constantelor şi , formula devine cos wt arctg cos arctg , w w - sau, dacă aplicăm proprietăţile funcţiilor trigonometrice, cos sin w w - Interpretare fizică Soluţia ne arată că amplitudinea oscilaţiilor scade exponenţial în timp Perioada mişcării este intervalul de timp care separă două treceri succesive printr-un extrem de acelaşi tip şi se determină impunând condiţia q& Cu ajutorul lui această condiţie se explicitează sub forma sinwt Rezultă astfel expresia perioadei oscilaţiilor amortizate vezi şi T g p p w g = × - l Ecuaţii diferenţiale ordinare de ordinul n Observăm că formula de mai sus se reduce la cazul pentru Este interesant de sesizat că, în ciuda amplitudinii descrescătoare, micile oscilaţii continuă să fie izocrone Aplicaţia Configuraţia de echilibru unui fir perfect flexibil solicitat concomitent de greutatea proprie şi de pretensionare orizontală G Cosovici, S Comşa Problema fizică Considerăm cazul unui fir solicitat de propria greutate şi de forţa orizontală H= figura Firul are următoarele caracteristici: lungimea l aria secţiunii transversale respectiv greutatea unităţii de lungime q , N/m Pentru ne face imagine asupra dimensiunilor secţiunii transversale, admitem că aceasta ar fi circulară de diametru d Folosind datele de mai sus, obţinem prin calcul m, A d p p - = × × deci un raport lungime diametru d » × × l Această din urmă valoare evidenţiază faptul că firul este foarte subţire În asemenea circumstanţe, poate fi adoptată ipoteza flexibilităţii perfecte, admiţând că firul este capabil să preia doar solicitări de întindere, fără să posede rezistenţă la încovoiere sau forfecare Figura Fir solicitat concomitent de greutatea proprie şi de pretensionare orizontală O y H x l q Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie O doua ipoteză modelului pe care îl vom prezenta mai jos poate fi justificată evaluând raportul dintre forţa de pretensionare şi greutatea întregului fir: H q = » l × Valoarea obţinută ne determină să admitem că variaţiile tensiunilor din fir cauzate de greutatea proprie sunt neglijabile în comparaţie cu efectul pretensionării prin forţa H Coroborând această simplificare cu ipoteza flexibilităţii perfecte, concluzionăm că eforturile de tracţiune au aproximativ aceeaşi componentă orizontală pe toată lungimea figura Model matematic Pentru descrie deformarea firului vom utiliza un sistem de coordonate carteziene cărui origine este amplasată la capătul din stânga al firului, axele şi fiind orientate aşa cum se vede în figura Valoarea foarte redusă a greutăţii proprii în comparaţie cu forţa de pretensionare ne determină să adoptăm ipoteza deformaţiilor mici, admiţând că secţiunile transversale ale firului se deplasează numai pe axa şi distorsiunile lor sunt neglijabile În aceste condiţii, deplasările verticale ale particulelor aflate iniţial pe axa Ox pot fi considerate reprezentative pentru toate celelalte particule din secţiunile transversale corespondente Altfel spus, săgeata verticală este funcţie numai de coordonata x: f xÎ[ l] Determinarea configuraţiei de echilibru firului se reduce la găsi funcţia f Problema se rezolvă impunând condiţia ca suma componentelor verticale ale forţelor care acţionează asupra unui element liniar de fir să fie zero figura : q dx tg da tga Trebuie observat faptul că tangenta unghiului este derivata lui în raport cu variabila x Ecuaţii diferenţiale ordinare de ordinul n d tg tg d d f f x x x a + Prin înlocuirea expresiilor ale tangentelor în egalitatea se obţine q dx dx sau, după eliminarea parantezei şi simplificarea cu , q f H Soluţie Ecuaţia este ecuaţie diferenţială lineară şi neomogenă, de ordinul al doilea Soluţia sa generală se obţine prin integrare directă de două ori şi are expresia + ÎÂ x c H q f Relaţia defineşte familie de parabole parametrizată de constantele şi Pentru obţinerea unei soluţii unic determinate, ecuaţia trebuie cuplată cu două condiţii suplimentare În cazul firului, respectivele condiţii impun anularea săgeţilor la capete figura : f Aplicând constrângerile de mai sus soluţiei generale deducem sistemul de ecuaţii algebrice liniare q c c H = Prin rezolvarea acestuia obţinem q c c H = l Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie După înlocuirea expresiilor ale constantelor de integrare, formula devine [ ] q f x H = Figura Element de fir deformat sub acţiunea propriei greutăţi şi unei pretensionări orizontale Aceasta este ecuaţia unei parabole, simetrică faţă de dreapta verticală / Potrivit relaţiei săgeata maximă corespunde abscisei / / şi are valoarea , m q q f H H × × l l l Aplicaţia Încovoierea grinzilor drepte Ecuaţia fibrei medii deformate G Cosovici, S Comşa Problema fizică Vom determina configuraţia de echilibru unei grinzi drepte solicitate la încovoiere pură figura Problema va fi rezolvată admiţând următoarele ipoteze: · Grinda are secţiunea transversală constantă pe toată lungimea sa O f a +df a +da x y x dx T H T +dT H q Ecuaţii diferenţiale ordinare de ordinul n · Centrele de greutate ale secţiunilor transversale sunt repartizate pe dreaptă care coincide cu axa Ox sistemului de coordonate originea fiind la capătul din stânga al grinzii · Secţiunea transversală este simetrică atât faţă de axa Oy direcţia verticală pe schiţa din figura cât şi faţă de axa Oz perpendiculara pe planul figurii · Dimensiunile secţiunii transversale sunt mici în comparaţie cu lungimea grinzii · Momentul încovoietor este constant pe toată lungimea grinzii şi are direcţia axei Oz · Materialul grinzii are comportare liniar elastică descrisă de legea lui Hooke · Deformaţiile de ansamblu ale grinzii sunt mici altfel spus, deplasările particulelor sale au valori foarte reduse în comparaţie cu dimensiunile caracteristice ale secţiunii transversale · În configuraţia de echilibru grinzii încovoiate, secţiunile transversale iniţial plane şi perpendiculare pe axa Ox rămân plane şi devin perpendiculare pe curba definită de noile poziţii ale centrelor de greutate secţionale postulatul lui Bernoulli Figura Grindă solicitată la încovoiere pură Model matematic Vom efectua analiza deformaţiilor de încovoiere procedând la separarea imaginară unui element infinitezimal de grindă figura Acest element este delimitat în configuraţia iniţială prin secţiunile transversale AB şi CD, având abscisele respectiv dx Aplicarea momentului determină curbarea l O x y M M Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie grinzii Altfel spus, segmentele AC şi BD îşi pierd rectiliniaritatea, transformându-se în arcele Ç A respectiv Ç B Axa centrelor de greutate ale secţiunilor transversale urmează aceeaşi evoluţie, curbându-se la rândul său Totuşi, potrivit postulatului lui Bernoulli, noile secţiuni şi îşi conservă planeitatea, fiind de asemenea perpendiculare pe noua linie centrelor de greutate secţionale În aceste condiţii, se poate considera că fibrele se curbează luând forma unor arce de cerc concentrice Deformaţia de ansamblu grinzii fiind mică, secţiunile transversale îşi vor conserva simetria Drept consecinţă, linia centrelor de greutate RS în figura îşi va păstra caracterul de fibră medie, iar lungimea sa nu va suferi modificări: = j Ç dx RS În relaţia au fost introduse următoarele notaţii: raza de curbură a fibrei medii Ç R dj unghiul la centru subîntins de fibrele longitudinale ale elementului de grindă deformat Fie mn fibră longitudinală elementului de grindă în configuraţie rectilinie poziţionată faţă de RS prin ordonata figura După aplicarea momentului M, segmentul mn se transformă în arcul de cerc Ç m Întrucât secţiunile grinzii nu îşi schimbă dimensiunile, se poate considera că raza lui Ç m este y Lungimea arcului Ç m este deci Ç m Cu ajutorul lui şi deformaţia fibrei mn dx se va explicita după cum urmează: Ecuaţii diferenţiale ordinare de ordinul n r = r j j e = Ç y y mn m mn d d d Figura Deformaţia unui element infinitezimal de grindă solicitat la încovoiere pură În stadiul imediat următor, legea lui Hooke, s Ee modulul de elasticitate permite evaluarea tensiunilor longitudinale datorate încovoierii: Ey s r = x dx x y y y r dj M M A B C D A B C D R S R m n m O Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Înlocuind expresia lui în condiţia de echilibru mecanic scrisă pentru momente M d S = domeniul plan reprezentat de secţiunea transversală grinzii vom obţine , EI unde const y S = este momentul de inerţie secţional raportat la axa Oz Atât timp cât deformaţiile grinzii sunt mici, deplasările particulelor aflate pe fibra medie au componentă orizontală neglijabilă În asemenea condiţii, raza de curbură r depinde numai de săgeţile verticale ale acestor particule Fie f xÎ[ ,l] funcţia care defineşte fibra medie în configuraţia deformată Dacă se ţine cont de orientarea sistemului de coordonate din figura raza de curbură este exprimabilă sub forma Ipoteza micilor deformaţii impune ca rotaţiile secţiunilor transversale să fie foarte reduse Acest fapt permite neglijarea termenului în conducând la expresia aproximativă r Prin înlocuirea lui definit de în obţinem ecuaţia diferenţială a fibrei medii: Ecuaţii diferenţiale ordinare de ordinul n z M f EI Soluţie Soluţia generală lui rezultă prin două integrări succesive în raport cu variabila x: + Î[ ÎÂ x c EI M f x z Relaţia defineşte parabolă parametrizată de constantele şi Pentru determinarea acestora, ecuaţia trebuie cuplată cu două condiţii care impun anularea săgeţilor de capăt: f Aplicând constrângerile soluţiei generale deducem sistemul de ecuaţii z M c c EI = Prin rezolvarea acestuia se obţine z M c c EI = l După înlocuirea expresiilor ale constantelor de integrare şi , formula devine [ ] z M f x EI = Această ultimă relaţie defineşte soluţia particulară ecuaţiei diferenţiale care verifică şi condiţiile la limită Formula de mai sus descrie parabolă simetrică faţă de dreapta verticală / De fapt, potrivit relaţiei săgeata maximă corespunde tocmai punctului de abscisă / şi are expresia Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie z M f EI × l l Aplicaţia Flambajul elastic al unei grinzi articulate la capete şi supuse unei solicitări axiale G Cosovici, S Comşa Problema fizică În anumite condiţii de solicitare, unele sisteme mecanice pot avea mai multe configuraţii de echilibru Figura ilustrează această posibilitate pentru cazul unei grinzi articulate la capete şi supuse unei solicitări axiale F Atât timp cât nivelul forţei nu depăşeşte un nivel critic Fcr grinda rămâne dreaptă preluând încărcarea în regim de compresiune pură Atunci când atinge valoarea Fcr oricare din cele două configuraţii de echilibru reprezentate în figura devine posibilă Teoretic, orice perturbaţie este în măsură să determine trecerea bruscă grinzii la forma curbilinie Acest fenomen se numeşte flambaj Este uşor de remarcat faptul că, în noua configuraţie de echilibru, încărcarea va fi preluată atât în regim de compresiune, cât şi prin încovoiere Dacă grinda este capabilă să revină la forma iniţială după eliminarea forţei flambajul se numeşte elastic În practică, pot fi întâlnite şi situaţii de flambaj ireversibil, atunci când solicitările din material depăşesc limita de curgere Figura Flambajul elastic al unei grinzi articulate la capete şi supuse unei solicitări axiale În ceea ce urmează vor fi stabilite condiţiile de trecere grinzii din figura de la configuraţia dreaptă la configuraţia curbilinie, admiţând că flambajul este de tip elastic După cum s-a menţionat anterior, în configuraţia sa curbilinie, grinda preia x x f y F l O Ecuaţii diferenţiale ordinare de ordinul n încărcarea axială atât în regim de compresiune, cât şi în regim de încovoiere Plecând de la această observaţie, flambajul va fi analizat scriind ecuaţia fibrei medii deformate şi căutând condiţiile în care ea admite mai multe soluţii Model matematic Considerăm secţiune transversală oarecare grinzii care şi-a pierdut forma rectilinie Această secţiune este poziţionată prin abscisa figura În condiţii de flambaj, săgeata fibrei medii defineşte un braţ al forţei axiale F Rezultă astfel un moment încovoietor M fF Printr-o metodă perfect similară celei de la aplicaţia precedentă, se poate deduce următoarea ecuaţie diferenţială care defineşte configuraţia de flambaj fibrei medii: z F f f EI Pentru comoditatea calculelor, este convenabilă definirea parametrului z F EI w × Cu ajutorul acestei mărimi se va rescrie sub forma +w Soluţie Ecuaţia este ecuaţie diferenţială ordinară de ordinul II, lineară şi omogenă, cu coeficienţi constanţi Soluţia ei generală este vezi ] sinw cosw xÎ ,l unde şi sunt constante care se determină impunând satisfacerea două condiţii la limită care reflectă blocajul mecanic exercitat de articulaţiile de la capetele grinzii figura : f Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Interpretare fizică Aplicând constrângerile funcţiei obţinem două ecuaţii în necunoscutele şi : sinw cosw sau, echivalent, sinw Flambajul grinzii poate fi descris numai de soluţie netrivială Pentru ca sistemul să devină compatibil cu această constrângere este necesar ca w să verifice egalitatea sinw Din rezultă condiţia w kp, ÎN sau, dacă ţinem cont de expresia parametrului , , , ÎN k l k EI F z cr Formula defineşte un şir infinit de valori critice ale forţei F Fiecare dintre aceste nivele ale încărcării produce un mod de flambaj Sub aspect practic, cea mai importantă este încărcarea minimă Aceasta corespunde indicelui : min z cr cr EI F F = × l Analizând relaţia constatăm că valoarea lui cr ,min este cu atât mai redusă cu cât raportul este mai mic altfel spus, cu cât dimensiunile transversale ale grinzii sunt mai mici în comparaţie cu lungimea Ecuaţii diferenţiale ordinare de ordinul n Aplicaţia Starea de tensiuni din tubul cilindric cu pereţi groşi G Cosovici, S Comşa Model matematic Considerăm cazul unui tub cilindric rectiliniu de lungime infinită, cărui secţiune transversală este un inel circular de raze figura Atât pe interiorul tubului, cât şi pe exterior acţionează presiuni presiunea interioară, respectiv presiunea exterioară Simetria axială tubului precum şi încărcărilor recomandă utilizarea unui sistem de coordonate cilindrice Axa coordonatei este coliniară cu axa de simetrie a tubului Cât priveşte coordonatele şi acestea se măsoară într-o secţiune transversală figura Lungimea infinită tubului, împreună cu simetria axială geometriei sale şi încărcărilor determină dependenţa deplasărilor, deformaţiilor şi tensiunilor exclusiv de raza r: , , , r z r z r z u r r r r r q q q q q q e e s s = = = = = = În egalităţile de mai sus au fost utilizate următoarele notaţii: componentele câmpului deplasărilor, deformaţiile principale şi tensiunile principale, asociate direcţiilor radială, circumferenţială, respectiv axială Datorită lungimii infinite tubului, deplasările şi deformaţiile axiale sunt nule: În plus, repartiţia simetrică încărcărilor determină anularea deplasării circumferenţiale: u Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Singura componentă diferită de zero câmpului deplasărilor este deci ur Aceasta determină complet deformaţiile εr şi e , r r r u u r = Admitem că materialul tubului are comportare elastică descrisă de legea lui Hooke: , , , r r r z r E E E q q q q s ne n n s ne e n n s e n n = = = + + în care este modulul lui Young, iar este coeficientul lui Poisson Figura Secţiune transversală printr-un tub cu pereţi groşi solicitat de presiunile şi pe suprafeţele sale interioară, respectiv exterioară Cu ajutorul relaţiilor formulele devin r r p p r s q s Ecuaţii diferenţiale ordinare de ordinul n r r r r r r z r E u u r E u u r E u u r q s n n n s n n n s n n n  +  Observăm că, în urma ultimei transformări, tensiunile au fost exprimate ca dependenţe de funcţia r Figura Schema de solicitare unui element de material separat din peretele tubului Aceasta din urmă se determină impunând condiţia de echilibrare radială a încărcărilor preluate de un element infinit mic separat din peretele tubului figura : d d sin d r r q s Folosind în egalitatea de mai sus aproximarea x y O r dr q dq d rs q d d d d r r s q dr s dr s Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie d d sin q şi operând serie de simplificări, obţinem +s Cu ajutorul formulelor condiţia de echilibru devine r u u u r r Soluţie Ecuaţia diferenţială este de tip Euler Ea are soluţia generală de forma vezi r A u Ar r = în care şi sunt constante Pentru determinarea acestora se impun condiţiile la limită semnul minus ia în considerare efectul de compresiune al încărcărilor aplicate asupra pereţilor tubului vezi figura Dacă se apelează la formulele şi , constrângerile de mai sus pot fi explicitate sub forma E A A p r E A A p r n n n n n n   Prin rezolvarea sistemului se obţin constantele p r A E r r r A p E r n n n + = - = - Ecuaţii diferenţiale ordinare de ordinul n După înlocuirea expresiilor de mai sus ale lui şi în rezultă soluţia particulară r p p u r E r +n +n × - În final, din şi deducem distribuţiile radiale ale tensiunilor din pereţii tubului: r z r p r r r p r r r p r r r q q s s s s - - - + - = + - Interpretare fizică Este interesant de remarcat faptul că tensiunea axială are aceeaşi valoare pe toată grosimea peretelui În ceea ce priveşte componentele r s şi s acestea au distribuţie variabilă pe rază, dar suma lor este de asemenea constantă Aplicaţia M V Soare [] Problema fizică Peste un catarg OA, de înălţime este trecut un cablu BOC figura Prin întinderea cablului se introduce forţă de compresiune în catarg Se cere să se determine valoarea forţei critice pentru care aceasta îşi pierde stabilitatea formei Model matematic Fie unghiul de înclinare al cablului faţă de orizontală, în poziţia iniţială figura Să presupunem că, datorită producerii flambajului, capătul superior suferă deplasare laterală f Atunci unghiul de înclinare părţii din stânga cablului se reduce cu un unghi Da, în timp ce unghiul de înclinare părţii din dreapta creşte cu Da figura Dacă este efortul de întindere în cablu, din condiţia iniţială de echilibru rezultă Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie a = sin P N Ca urmare deformării ansamblului, se va dezvolta forţă orizontală H cos Da cos Da sin asin Da; cum Da este foarte mic în comparaţie cu sin Da Da , H sin aDa PDa Dacă este proiecţia punctului pe dreapta BO este punctul în care ajunge O prin flambaj atunci din triunghiul ODO rezultă figura OD BODa sin a; deoarece BO sin se obţine = a a Da sin sin sin l f l f , Figura Schema geometrică catargului şi cablului; deplasarea laterală f; deplasarea capătului superior; schema statică catargului aşa încât Ecuaţii diferenţiale ordinare de ordinul n = sin a l f H Avem aşadar de determinat forţa critică unei console OA, încărcată la capătul liber cu forţele şi figura Alegând originea x-ilor la capătul superior al barei, ataşată acesteia, într-o secţiune curentă momentul încovoietor se scrie sin a l f M Pw Hx Ecuaţia diferenţială fibrei medii deformate barei va fi a sin d d l f w x EI P x w , unde EI este rigiditatea la încovoiere, sau + a sin d d l f w x x w cu notaţia uzuală EI P Soluţie Soluţia generală ecuaţiei lineare de ordinul II, cu coeficienţi constanţi se scrie vezi = sin cosb sin a l f w iar rotaţia secţiunii transversale este dată de = cosb sin sin a d d l f A x x w Condiţiile la capete , x=l conduc la următorul sistem de ecuaţii lineare şi omogene în şi f: Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie cos sin sin sin cos sin , , b = b = = l f A l A f B Ţinând seama că deducem  =   b a b a sin cos sin cos f A l l l Acest sistem are soluţii nenule corespunzând poziţiei stabile de echilibru dacă sin cos sin cos det =  b a b a l l l ; se obţine astfel următoarea ecuaţie caracteristică = a b b cot tan l l , ale cărei soluţii pot fi determinate numai numeric De exemplu, pentru se obţine bl , deci forţa critică este dată de l EI l EI Pcr = Aplicaţia Flambajul barelor Probleme de mecanica construcţiilor M V Soare , Flambajul barei dublu articulate Problema fizică bară zveltă este comprimată axial la capete de două forţe La atingerea unei valori critice forţelor Pcr bara părăseşte forma rectilinie de echilibru Ştiind că bara este dublu articulată la capete, se cere să se determine forţa critică de flambaj primele două valori şi forma fibrei medii deformate barei Ecuaţii diferenţiale ordinare de ordinul n Model matematic Pentru rezolvarea problemei, se consideră momentul în care bara părăsit forma rectilinie de echilibru şi căpătat formă curbilinie, foarte apropiată de poziţia iniţială În această situaţie, într-o secţiune curentă momentul încovoietor va fi Pw, unde este deplasarea, astfel încât ecuaţia diferenţială a fibrei medii deformate barei se scrie vezi şi aplicaţia EI Pw EI M x w d d adică d d + = EI P x w în care EI este rigiditatea la încovoiere minimă secţiunii transversale barei Pentru simplitate, folosim din nou notaţia şi ecuaţia devine d d + = x w Alegând originea x-ilor la capătul superior, condiţiile la limită sunt w unde este lungimea barei figura Soluţie Modelul reprezintă, din punct de vedere matematic, o problemă bilocală omogenă, care admite întotdeauna soluţia identic nulă Trebuie deci determinate valori ale lui denumite valori proprii, astfel încât problema să admită cel puţin încă soluţie, diferită de cea identic nulă, denumită funcţie proprie Asemenea probleme poartă numele de probleme Sturm-Liouville Căutând soluţii de forma elx vezi ajungem la ecuaţia caracteristică l cu rădăcinile ib Soluţia generală ecuaţiei se poate scrie sub forma Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie w Asin bx cosbx A şi fiind constante de integrare Condiţiile la limită implică şi Asin bl Cum axa barei este rectilinie, deci situaţia premergătoare flambajului iar b bara este neîncărcată rezultă sin bl cu rădăcinile bl np, K Revenind la notaţia se obţine familia de valori proprii , K = = l EI Pcr EI n şi ecuaţia fibrei medii deformate, care se deduce din ţinând seama de valorile găsite pentru şi b sin n l n x w Se observă că factorul reprezentând amplitudinea fibrei medii deformate, rămâne nedeterminat Figura Flambajul barei dublu articulate Ecuaţii diferenţiale ordinare de ordinul n Aceasta se explică prin faptul că, la rezolvarea problemei, fost utilizată ecuaţia diferenţială aproximativă fibrei medii deformate; ecuaţia reprezintă o sinusoidă cu semiunda n Din punct de vedere practic, interesează valoarea minimă forţei critice pentru n Aceasta mai este numită şi forţa critică Euler-iană l EI Pcr PE = Fibra medie deformată este sinusoidă de semiundă şi este dată de max l x w w = unde wmax corespunde mijlocului deschiderii Pentru valori mai mari ale lui de exemplu, aceasta ar corespunde cazului când există la mijlocul deschiderii, un reazem simplu Forţa critică corespunzătoare este cr Pcr l EI P , Flambajul barei articulate-încastrate Problema fizică Să se determine forţa critică de flambaj unei bare articulate la un capăt şi încastrată perfect la celălalt capăt Model matematic Datorită încastrării, în articulaţie apare şi reacţiune H, normală pe axa barei, jucând rolul unui parametru nedeterminat Momentul încovoietor într-o secţiune curentă evaluat pe forma deformată, este Pw Hx astfel încât ecuaţia diferenţială problemei este x EI H w EI P x w d d unde este forţa de compresiune şi EI rigiditatea la încovoiere Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Figura Flambajul barei articulate-încastrate Cu notaţia ecuaţia devine x EI H w x w d d Condiţiile la limită se scriu w , x=l Soluţie Modelul de mai sus reprezintă tot problemă Sturm-Liouville Soluţia generală ecuaţiei omogene asociate lui este tot soluţie particulară rezultă sub formă de monom Soluţia generală ecuaţiei neomogene şi derivata sa se scriu sub forma cos sin d d sin cos , A x P H x w x x P H w = b = b Ecuaţii diferenţiale ordinare de ordinul n Prima condiţie la limită implică Celelalte două condiţii conduc la un sistem algebric linear, pe care trebuie să-l satisfacă şi A: - sin Acosbl P H l l P H Acest sistem este omogen, deci are soluţii nenule doar dacă determinantul asociat se anulează cos sin det  b b b l l l l Calculând determinantul, obţinem ecuaţia transcendentă tan bl bl Rădăcina minimă acestei ecuaţii conform tabelului b l conduce la forţa critică minimă l EI Pcr @ Flambajul barei dublu încastrate Problema fizică Să se determine forţa critică de flambaj unei bare articulată la un capăt şi încastrată perfect la celălalt capăt Model matematic Ţinând seama şi de cazul precedent, ecuaţia diferenţială a problemei este EI M x EI H w EI P x w d d unde este forţa de compresiune, şi reacţiunea din origine, normală pe axa barei, respectiv momentul de încastrare sunt parametri nedeterminaţi, EI este rigiditatea Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie la încovoiere şi este săgeata necunoscută figura Condiţiile bilocale se scriu în acest caz astfel w , x= x=l unde este lungimea barei Figura Flambajul barei dublu încastrate Soluţie Modelul matematic se prezintă sub forma unei probleme Sturm-Liouville Mai întâi găsim soluţia ecuaţiei lineare şi neomogene vezi aceasta şi derivata ei sunt cos sin d d sin cos , A x P H x w A x P M x P H w = b = b Condiţiile bilocale conduc la sistemul algebric linear şi omogen, scris sub forma matriceală Ecuaţii diferenţiale ordinare de ordinul n  =   b b b - cos sin sin cos M P H P B A l l l l Pentru ca sistemul să admită soluţii nenule trebuie ca determinantul asociat lui să se anuleze Se obţine astfel ecuaţia caracteristică tan sin sin cos - bl l l ale cărei rădăcini sunt bl np K Rădăcina cea mai mică conduce la forţa critică l EI l EI Pcr = Rădăcina minimă corespunzătoare celui de al doilea factor este mai mare decât b Bara dublu articulată încărcată cu forţe axiale şi cu forţă transversală Problema fizică bară dreaptă, dublu articulată, de lungime este solicitată axial prin două forţe de compresiune şi încărcată transversal cu forţă concentrată , aplicată la mijlocul deschiderii figura Se cere să se determine săgeţile w Model matematic Pentru secţiune transversală de abscisă putem scrie momentul de încovoiere M Pw Fx Problema fizică este deci modelată de ecuaţia diferenţială lineară de ordinul II cu coeficienţi constanţi Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie EI l P x x EI F w x w  , d d , căreia se asociază condiţiile la capete x=l ultima este condiţie de simetrie Modelul matematic este reprezentat de problemă bilocală lineară Figura Bara dublu articulată încărcată cu forţele axiale şi cu forţă transversală F Soluţie Soluţia generală ecuaţiei de mai sus este vezi x x P F w sin cosb , şi deci A x P F x w cosb sin b d d Aplicând condiţiile la capete, rezultă , cos = b b = B l P F A , astfel încât b - b b = cos sin l x x P F w Momentul de încovoiere are expresia Ecuaţii diferenţiale ordinare de ordinul n cos sin l F x M b b b = Pentru P® PE EI forţa lui Euler avem cos bl deci şi tind la infinit, independent de intensitatea forţei transversale instabilitate prin divergenţă Aplicaţia M V Soare , Problema fizică Să se studieze flambajul lateral al unei grinzi cu secţiune zveltă, supusă la încovoiere Model matematic Grinzile supuse la încovoiere având secţiune zveltă îşi pot pierde forma plană de echilibru atunci când momentul încovoietor depăşeşte anumită valoare critică figura Grinda îşi pierde stabilitatea în zona comprimată; axa grinzii se curbează în planul ei de rigiditate minimă, iar diferitele secţiuni transversale ale grinzii se rotesc în jurul axei Acest fenomen de pierdere stabilităţii formei de schilibru unei grinzi supuse la încovoiere este denumit flambaj lateral sau flambaj din încovoiere Studiul flambajului lateral conduce la ecuaţia diferenţială ordinară lineară şi omogenă d d q = EI zGIt M x în care · reprezintă rotirea de torsiune secţiunii transversale · EI rigiditatea la încovoiere minimă grinzii în planul z, · GIt rigiditatea la torsiune secţiunii grinzii, · momentul încovoietor în planul y Sunt considerate secţiunile simple fără tălpi Introducând notaţia Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie EI zGIt M ecuaţia devine d d q = x analogă celei corespunzătoare flambajului barei dublu articulate de mai sus Acestei ecuaţii se asociază condiţiile bilocale q obţinându-se astfel problemă Sturm-Liouville Figura Flambajul lateral al unei grinzi cu secţiune simplă Soluţie Soluţia generală ecuaţiei este q Asin bx cosbx Aplicând condiţiile la capete, se obţine valoarea proprie minimă aşa încât cr EI zGIt l M = Ecuaţii diferenţiale ordinare de ordinul n Aplicaţia M V Soare , Problema fizică bară de oţel, în consolă, încărcată cu forţa este rezemată lateral cu un resort având constanta elastică c Se cere să se determine forţa critică de flambaj Pcr Model matematic Expresia momentului încovoietor într-o secţiune transversală curentă de abscisă este dată de figura M cf unde este forţa axială, deflexia la capătul cu reazem elastic constanta elastică este c şi este lungimea barei Figura Bara încărcată axial cu forţa încastrată şi rezemată la cele două capete Cu notaţia rezultă ecuaţia diferenţială fibrei medii deformate barei  x P cf w f x w d d la care se adaugă condiţiile la limită x= Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie deci din nou problemă de funcţii şi valori proprii Soluţie Soluţia generală ecuaţiei diferenţiale este x P cf P cl f sin cosb Condiţiile iniţiale de la marginea inferioară conduc la EI cl B EI cl A f , b = b = astfel încât pentru fibra medie rezultă expresia  b b - b + b = x EI c x EI cl w sin cos Condiţia duce la ecuaţia caracteristică sin cos b b - b + b - l EI c l EI cl care se mai poate scrie bl bl tan bl unde cl EI k Pentru bară de secţiune circulară de diametru şi datele numerice E × daN cm daN cm, cm, obţinem = × × p× k Ecuaţii diferenţiale ordinare de ordinul n Tabelul Valorile lui în funcţie de βl bl bl k π/ +¥ Pentru această valoare constatei rădăcina minimă ecuaţiei este bl Forţa critică devine l EI l EI Pcr = , unde fost pusă în evidenţă şi lungimea de flambaj Figura Graficul funcţiei bl Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie În cazul altor date numerice, în tabelul au fost date valorile lui pentru diferite valori bl vezi relaţia : b - b b = b = b l l l l l l k l tan tan Graficul funcţiei bl este redat în figura Atât tabelul cât şi figura pot servi la determinarea rădăcinii bl pentru un dat Aplicaţia M V Soare , Problema fizică Să se studieze flambajul unei grinzi drepte într-un mediu elastic, în cazul general de rezemare, în ipotezele lui Winkler kw, compresiunea este proporţională cu deplasarea const fiind coeficientul reacţiei solului , Model matematic Ecuaţia diferenţială care guvernează deformarea grinzii este d d + = j w w unde parametrul depinde de elasticitatea mediului Soluţie Ecuaţia este lineară şi cu coeficienţi constanţi vezi Căutând soluţii sub forma exponenţială ebx se obţine ecuaţia caracteristică l , cu rădăcinile complex conjugate Soluţia generală poate fi exprimată în una din următoarele forme: cosh sin sinh sin , cosh cos sinh cos A x w x + b = b w e-bx Acosbx sin bx ebx cosbx Dsin bx în care respectiv B,C,D, sunt constante de integrare Ecuaţii diferenţiale ordinare de ordinul n Plecând de la formularea vom schimba constantele de integrare, introducând altele noi, care au semnificaţie fizică parametri iniţiali ,j ,M ,T , reprezentând săgeata, rotirea, momentul încovoietor şi forţa tăietoare la capătul din stânga al grinzii ales ca origine x-ilor Figura Graficele funcţiilor bx i Dacă se notează funcţiile sinh cos cosh sin , sinh sin , sinh cos cosh sin , cosh cos , f x f x f x f x b b b b b b b b expresiile săgeţii, rotirii, momentului încovoietor şi forţei tăietoare se scriu: Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie d d , d d , d d , f x k f x kw x w T EI f x T f x k f x kw x w M EI f x k T f x k M w x x w f x k T f x k M w x b b b b j b = = b b b + b b j b = = j b b b b = j Observaţie Funcţiile fi bx definite de au largă răspândire în mecanica construcţiilor Graficele lor sunt redate în figura Aplicaţia M V Soare , Problema fizică grindă de lungime foarte mare teoretic infinită reazemă pe un mediu elastic şi este încărcată cu forţa transversală P se cere să se determine expresia săgeţii rotirii momentului încovoietor şi forţei tăietoare într-o secţiune curentă, apoi să se reprezinte grafic Model matematic Grinda fiind infinită, originea axelor poate fi aleasă în orice punct; va fi avantajos să alegem în punctul de aplicare forţei concentrate, pentru a profita de condiţiile de simetrie şi antisimetrie figura La aplicaţia precedentă aplicaţia s-a găsit expresia generală săgeţii w e-bx Acosbx sin bx ebx cosbx Dsin bx unde este constantă dată de / EI fiind coeficientul de tasare al mediului elastic şi EI rigiditatea la încovoiere grinzii Ecuaţii diferenţiale ordinare de ordinul n Figura Grinda pe mediu elastic, de lungime infinită: diagrama săgeţilor; diagrama forţelor tăietoare în vecinătatea originii Soluţie În origine vom avea dw/ dx iar forţa tăietoare are un salt; la dreapta, Tr figura La infinit, unde efectul forţei exterioare nu se mai resimte, putem exprima că orice triplet de mărimi şi tinde la zero Com factorul ebx creşte nedefinit când x®¥, va trebui să luăm Rămân astfel expresiile [ ] cos sin ] d d cos sin , d d e cos sin , d d e cos sin , EI x x w T EI EI x x w M EI B x x w w x x x x x = b = b j b = b -b -b -b -b Introducând condiţiile la limită precedente, rezultă Pb Cu aceasta, expresiile devin cos pentru e cos sin , e cos sin sin , b b = = b -b b -b b x x P x P M x k P x x k P w x x x x Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Figura Graficele funcţiilor w,j,M,T Se observă că în formulele apar patru funcţii de argument bx numite funcţii exponenţiale amortizate cos sin e cos sin , e sin , e cos , x x x x x x x x x x x x y b y b y b y b -b -b -b -b Observaţie Funcţiile yi bx au şi ele largă răspândire în mecanica construcţiilor Graficele lor sunt redate în figura Ecuaţii diferenţiale ordinare de ordinul n Aplicaţia M V Soare , Problema fizică particulă electrizată intră intr-un câmp electromagnetic de intensitate şi inducţie B Se cere să se determine traiectoria particulei Model matematic Componentele celor două forţe faţă de un triedru ortpgonal de referinţă Oxyz sunt reprezentate în figura Forţa rezultantă este qE qv B, unde este sarcina electrică, iar al doilea termen este forţa Lorenz Avem încă i j i k v B B = vx vz x x cu v vxi vy vzk x&i y&j z&k pentru studiul mişcării, introducem legea lui Newton F ma m&r& unde este masa particulei Soluţie Proiectând pe cele trei axe de coordonate se obţin ecuaţiile de echilibru dinamic m&x& qv yB m &y& qEy qvxB m&z& qEz Ultima ecuaţie poate fi considerată separat Integrând dată, obţinem d m qE t m qE z& ∫ Constanta se determină din condiţia iniţială z& vz Rezultă C vz şi Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie z z v m qE z& O nouă integrare dă t C m qE z z = Figura Particula electrizată într-un câmp electromagnetic Condiţia implică astfel încât t t m qE z z z = reprezintă mişcare uniform accelerată în lungul axei Oz, de acceleraţie az qEz m Pentru celelalte două axe, ţinând seama de ecuaţiile şi se scriu m&x& qBy& &y& qB qBx& Eliminând funcţia din prima ecuaţie şi introducând-o în cea de doua, deducem m q BE x m q B x &&& & Ecuaţii diferenţiale ordinare de ordinul n deci ecuaţie diferenţială ordinară de ordinul III, lineară şi neomogenă, cu coeficienţi constanţi vezi Raportul reprezintă sarcina electrică unităţii de masă Cu notaţiile = w m q B , R m q BEy ecuaţia devine &x&& Cum termenul liber este constantă, găsim uşor soluţie particulară acestei ecuaţii xp Eyt B Pentru ecuaţia omogenă asociată &x&& se caută soluţia sub forma exponenţială elt şi rezultă ecuaţia caracteristică l , cu trei rădăcini iw Soluţia generală ecuaţiei omogene asociate este deci xh coswt sin wt În final, soluţia generală ecuaţiei neomogene sau, echivalent, este t B E x cosw sin Din prima ecuaţie rezultă anume const cos sin const = wt wt B E qB m x qB m y y & adică Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie = sin wt coswt const qB m y Cele patru constante de integrare se determină din condiţiile iniţiale , y x& vx y& Se găsesc pentru ele următoarele expresii w = m qBv D , w = m qBv D , w - w = B v E D y , qB mE qB mv D , astfel încât cos sin cos sin sin , t v t qB mE y t v t t B E t m qBv x y y y x w w = + w w w + w - + w = Ţinând seama de notaţia deplasările şi se scriu sin cos sin cos sin , t x B E v t qB m y t t B E v t qB m x y y y x y &  = w  = w Am obţinut astfel ecuaţiile parametrice ale proiecţiei traiectoriei particulei electrizate pe planul xOy; curba obţinută se numeşte trohoidă CAPITOLUL SISTEME DE ECUAŢII DIFERENŢIALE ORDINARE Forma generală unui sistem de EDO cu funcţii necunoscute este , , , = = = n n n n n n F y F y F y unde Fj sunt definiţi pe un acelaşi domeniu dimensional şi sunt suficient de regulaţi Dacă ipotezele teoremei funcţiilor implicite sunt satisfăcute vezi cursul de Analiză Matematică, partea atunci pot fi explicitaţi din se obţine astfel forma canonică/normală unui sistem de EDO de ordinul I: , , , , n n n n y y y y y y = = = Sistemele de forma pot fi scrise şi în formă compactă Într-adevăr, utilizând notaţiile , , , , d d ,  =  =  = y y y f y y y f x f x f x x y y y x y y y n n M M sistemul poate fi scris în forma vectorială Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie f y , d d x x = Să mai observăm că orice ecuaţie diferenţială ordinară de ordinul n y y, , n- poate fi scrisă sub forma unui sistem de EDO de ordinul I Cu notaţiile , , = n- y yn ecuaţia devine , , , , n n n n y y y y y y y y = = = = - care este un sistem de EDO de ordinul în funcţii necunoscute Se poate demonstra că, invers, un sistem de EDO de ordinul de formă canonică poate fi redus la EDO de ordinul în anumite condiţii de regularitate SISTEME DE EDO DE ORDINUL LINEARE Forma canonică unui sistem neomogen de EDO lineare de ordinul în funcţii necunoscute, este , , , y x y x y x n nn n n n n n + + + + + + unde fi aij sunt funcţii considerate de clasă [a,b]ÎÂ iar y sunt funcţiile necunoscute Sisteme de ecuaţii diferenţiale ordinare Lucrăm într-un cadru clasic, deci vom căuta soluţii cel puţin netede, adică de clasă C Dacă fi se anulează identic pe obţinem sistemul omogen asociat , , n nn n n n n n y y y y y y + + + Pentru simplificarea scrierii şi, de asemenea, pentru pune în evidenţă unele aspecte utile în aplicaţii, vom introduce următoarele funcţii vectoriale:  =  =  = f x f x f x x y x y x y x x y x y x y x x n n M M , d d , f y y precum şi matricea asociată lui sistemului sau :  = a x a x a x x n nn n n K K K K K A Cu aceste notaţii, sistemul neomogen se scrie sub forma vectorială A f y = + dx d iar sistemul omogen asociat ia forma A y y x x = d d Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Un prim rezultat important, care este consecinţă directă linearităţii, afirmă că orice combinaţie liniară de soluţii ale sistemului omogen este de asemenea o soluţie sistemului Teorema Fie matrice , funcţie continuă de x Dacă şi y sunt două soluţii ale sistemului omogen atunci C +C ÎÂ, este de asemenea soluţie sistemului Al doilea rezultat afirmă că soluţia generală sistemului neomogen se obţine ca suma dintre soluţie particulară lui şi soluţia generală sistemului omogen Teorema Fie matrice şi funcţii continue de x Notăm cu o soluţie particulară sistemului neomogen Atunci este soluţia generală a sistemului dacă şi numai dacă yo unde yo este soluţia sistemului omogen Demonstraţie Observăm că ' ' x x x x p o p o A f y y = + + = Reciproc, dacă yo este soluţia sistemului neomogen, atunci , ' ' x x x x p p p A y y f = - = adică este soluţia sistemului omogen Sisteme de ecuaţii diferenţiale ordinare SISTEME DE EDO DE ORDINUL LINEARE, CU COEFICIENŢI CONSTANŢI Teoremele generale şi nu descriu modalitatea de determinare explicită a soluţiilor sistemelor lineare, omogene sau neomogene Vom vedea, în cele ce urmează, că putem obţine soluţii explicite dacă matricea este constantă nu depinde de Forma generală unui sistem linear de ordinul cu coeficienţi constanţi este        + + + + + + , , , y x y x y x n nn n n n n n unde ÎC , ÍÂ n, aij ÎÂ,i Cu notaţiile sistemul neomogen se scrie vectorial astfel: Ay f y + dx d unde matricea constantă este dată de = n nn n n a a a a a a A La fel ca mai sus, pentru obţinem sistemul linear omogen asociat Ay y = dx d Pentru rezolvarea sistemului căutăm soluţii de forma cn n e , = ÎÂ a C Derivând în raport cu obţinem Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie C y x x = aea d d Înlocuind în obţinem aeaxC eaxAC Simplificăm cu eax şi notăm cu matricea unitate: Î n´n E , M Obţinem relaţia aE unde este vectorul identic nul, cu componente Dacă sistemul admite soluţii de forma atunci  este valoare proprie matricei sistemului, iar  este vectorul propriu corespunzător Condiţia reprezintă, de fapt, un sistem algebric linear şi omogen de n ecuaţii, cu necunoscute , , cn scris matriceal Problema rezolvării sistemului diferenţial omogen se reduce astfel la determinarea vectorilor şi valorilor proprii ale matricei asociate Pentru ca să admită soluţii care să nu fie identic nule, trebuie ca det aE adică, mai explicit, Sisteme de ecuaţii diferenţiale ordinare n n n nn a a a a a a - a - a = Condiţia sau, echivalent reprezintă ecuaţie algebrică polinomială, care se numeşte ecuaţie caracteristică De soluţiile ei depinde forma soluţiei sistemului de EDO Pentru exemplificare, vom considera cazul :        = + = + = + d d , d d , d d a y x y a y x y a y x y Ecuaţia caracteristică este a a a a a a - a - = - a şi are trei rădăcini ,a Distingem trei cazuri: A Cazul rădăcinilor reale şi distincte Fie ,a ,a ÎÂ distincte Atunci este valabilă schema e x a a ¯ ¯ a a Fie ,C ,C vectorii proprii corespunzători lui ,a ,a Rezultă, conform celor spuse anterior, că } ea xC ,ea xC ,ea xC formează un sistem fundamental de soluţii Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Soluţia generală sistemului este deci x y × × × unde ,k ,k sunt constante arbitrare B Cazul rădăcinilor complexe Dacă ecuaţia admite soluţie complexă, întrucât ea are coeficienţi reali, va admite şi conjugata ca soluţie Fie deci ib, ib, rădăcini complexe conjugate, şi ÎÂ Fie iD vectorul propriu asociat lui ib Atunci iD este vectorul propriu pentru ib vezi cursul de Algebră În loc să considerăm funcţiile complexe care rezultă direct, vom lua combinaţiile lineare ale acestora Deci C iD iD xC i i + a+ a- a ¯ a ¯ a b ¯ a b Atunci semisuma primelor două soluţii este de asemenea soluţie a sistemului, deoarece acesta este linear Avem [ ] , i e e e e e e e i e e e i e sin i i cos i i i i i i  = =  + = = b b b b b b a b b a a+ a- b x x x x x x x x x x x x C D C D y D de unde Sisteme de ecuaţii diferenţiale ordinare y eax Ccosbx Dsinbx Aceasta este, de fapt, partea reală lui C+ iD a+ib Obţinem altă soluţie reală scăzând primele două soluţii din şi împărţind această diferenţă la i: [ ] , e e i e e e i e e i i e e e i e i cos i i sin i i i i i i  = + =  + = = b b b b b b a b b a a+ a- b x x x x x x x x x x x x C D C D y D deci y eax Csinbx Dcosbx care este partea imaginară lui C+ iD a+ib În final, soluţia generală sistemului se scrie în acest caz e e cos sin sin cos C y D x x x k k x a a a + = + C Cazul rădăcinilor multiple Dacă ecuaţia admite rădăcină dublă cea de treia fiind diferită de ea, adică atunci schema sistemului fundamental eaxC ea xC ¯ ¯ a a este, deocamdată, incompletă Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie În diagrama de mai sus, şi sunt vectorii proprii corespunzători lui , respectiv Mai este necesar un element pentru sistemul fundamental Ca şi în cazul similar al EDO lineare căror ecuaţie caracteristică admite rădăcini multiple, căutăm soluţii de forma y xCeax Deax Înlocuim în şi determinăm pe din sistemul algebric ce se obţine după simplificarea cu eax Soluţia generală sistemului este k x y + Exemple Pentru cazul Să se rezolve sistemul        = = d d , d d y x z z x y Soluţie Este un sistem de ecuaţii diferenţiale linear şi omogen Matricea asociată sistemului este = A Notând cu = z y y vectorul funcţiilor necunoscute şi cu = z y y vectorul derivatelor lor, scriem sistemul sub formă matriceală y Ay Căutând soluţii de forma eaxC, ajungem la ecuaţia caracteristică Sisteme de ecuaţii diferenţiale ordinare det = - a - a A aE ale cărei soluţii sunt reale şi distincte: ,a Vectorii proprii sunt - = = , C Deci x x - = = e e , y formează bază sistem fundamental de soluţii pentru sistemul Soluţia sa generală este x - + = e e y sau, pe componente, ÎÂ    = = + - - e e e e k k z k y k x x x x , arbitrare Pentru cazul Să se rezolve sistemul        = = d d , d d y x z z x y Soluţie Matricea asociată sistemului este - = A Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie şi, utilizând aceleaşi notaţii ca în exerciţiul precedent, obţinem forma matriceală a sistemului: y Ay Ecuaţia caracteristică det A E -a - = - a are rădăcinile pur imaginare: Vectorii proprii corespunzători sunt complex conjugaţi - = = i , i C Rezultă sistemul fundamental ix ix i e , i - = Y Calculăm semisuma soluţiilor şi semidiferenţa împărţită la i Obţinem cos sin e e i e e i Im , sin cos i e e e Re i i i i i i i i = + - = = - = - + = = - - - - x x x x x x x x x x x x Y Y y Y Y Y y Y Soluţia generală sistemului este deci + - = x x k x x k cos sin sin cos y sau, pe componente Sisteme de ecuaţii diferenţiale ordinare cos sin sin cos y x k k z x =   ÎÂ =  , arbitrare Pentru cazul Să se rezolve sistemul        = = + d d , d d az x z ay z x y Soluţie Folosind matricea asociată = a a Α şi aceleaşi notaţii pentru funcţiile necunoscute şi derivatele lor, scriem sistemul sub forma matriceală y Αy Ecuaţia caracteristică este det a A a a - a - = - a având rădăcina dublă Vectorul propriu corespunzător lui este = C de unde obţinem primă soluţie sistemului fundamental căutat eax y Căutăm pe de forma Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie ax ax d d e e + y Înlocuim în sistem: ax ax ax ax ax d d a a x a a d d ax e e e e e + = + + Simplificăm cu eax şi rezultă + = + ad ad d ad ad de unde obţinem , fiind arbitrar Putem lua Astfel, cea de doua soluţie a sistemului fundamental este ax eax e + y sau x ax e y Soluţia generală sistemului este ax ax k e e + y sau, pe componente ÎÂ    = = + e e k k z k y x ax ax , arbitrare Observaţii Constantele arbitrare din soluţia generală sistemului se determină din condiţii suplimentare De multe ori, acestea sunt de tip Cauchy, sau iniţiale: Sisteme de ecuaţii diferenţiale ordinare y , , yn yn sau, altfel scris, unde T y , yn Pentru sistemul neomogen observăm, ca şi în cazul EDO lineare, că soluţia generală se scrie ca y yomogen unde este soluţie particulară lui iar yomogen soluţia generală sistemului omogen asociat Această observaţie este valabilă şi pentru sistemul general Soluţia particulară se caută fie utilizând metoda variaţiei constantelor, valabilă şi în cazul general, fie, ca şi în cazul EDO lineare cu coeficienţi constanţi, de forma termenului liber Exemplu Să se rezolve sistemul        = = + d d e , d d y x z z x y x Soluţie Sistemul omogen asociat este chiar cărui soluţie generală am găsit-o anterior Deci k x - + = e e yomogen Căutăm soluţie particulară sistemului neomogen de forma x b a Y Introducem în sistem şi obţinem Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie + = e e e x x x b a b a Simplificând cu , deducem sistemul algebric în b , , b a a b = = + deci , Rezultă e x Y astfel încât soluţia generală sistemului neomogen este k x e e e + - + y EXPRIMAREA SOLUŢIEI UNUI SISTEM DE EDO LINEARE FOLOSIND EXPONENŢIALA DE MATRICE Fie acum Mn spaţiul Banach al matricelor de numere reale înzestrat cu norma ΣΣ = = = n i n j aij A , pentru AÎMn aij jÎ{ ,n} Definiţia Fie Suma seriei convergente Σ ¥ = ! k k k A se numeşte exponenţiala matricei şi se notează cu eA Sisteme de ecuaţii diferenţiale ordinare Dacă matrice are pe diagonală valorile n l ,l , , iar restul elementelor sunt egale cu zero, vom scrie diag ,l , ,ln Proprietăţile exponenţialei sunt: a n e , elI elI pentru orice lÎÂ; b dacă diag ,l , ,ln atunci eA diag el ,el , ,eln ; c dacă AB BA atunci eAeB eBeA eA+B ; d pentru orice matrice CÎMn cu determinantul nenul avem A ACC - e ; e tÎÂ t et et , d d Teorema Fie matrice constantă Sistemul , d d Ay y y = x admite unică soluţie dată de expresia y exAy Demonstraţie Observăm că n e şi y' AexAy Ay Unicitatea rezultă din teorema Cauchy-Lipschitz Sistemul Ay y = dx d admite ca soluţie generală y exAC, unde este un vector constant de componente , C Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Teorema Fie un interval din care conţine punctul matrice n constantă şi funcţie continuă pe J Sistemul de EDO neomogen d d Ay y y = x admite unică soluţie dată de formula du x x ∫ y A Demonstraţie Sistemul de EDO omogen , d d Ay y y = x admite soluţia unică xA Soluţia sistemului neomogen se caută sub forma y =G Avem ' ' y Ay AGc f y G'c Gc' AGc Gc' = + = + Urmează că satisface relaţia c' G- care admite soluţia unică x du x∫- + = f G y c În acest fel, x du x x Ay ∫e Af În teoremele precedente nu se calculeaza explicit exponenţiala, acest lucru fiind dificil, mai ales în cazul sistemelor de dimensiuni mari Dacă matricea este diagonalizabilă, calculul exponenţialei se face simplu Sisteme de ecuaţii diferenţiale ordinare Teorema Fie matrice n´ constantă, diagonalizabilă Atunci exA PexDP- unde este matricea diagonală ale cărei elemente sunt valorile proprii ale lui şi P este matricea ale cărei coloane sunt vectori proprii ai lui A Demonstraţie fiind diagonalizabilă, P- AP şi An PDnP- cu Dn matrice diagonală Se deduce că exD este matrice diagonală cu elementele exl ,exl , ,exln şi exA PexDP- Exemplu Să se rezolve sistemul , d d Ay y = x unde A= - - - Soluţie este diagonalizabilă şi valorile proprii sunt: = cu vectorii proprii , , = = v v Atunci - - - = - = = , , P , e e e e e e e = = - - - - x x x x x x xA xDP Soluţia generală sistemului este Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie e e e e x x x x x x y C y C y C = + = + = + - - - Dacă matricea nu este diagonalizabilă, pentru calcula exA se poate folosi următoarea teoremă Teorema Fie matrice constantă Atunci are unică descompunere unde este diagonalizabilă, este nilpotentă adică pentru k suficient de mare Nk iar şi comută Deoarece şi comută, deducem că exA exSexN fiind diagonalizabilă, exS se calculează ca mai sus Pentru exN avem: ! e - - = k k n x k x N xN SISTEME DE ORDINUL NELINEARE SISTEME SIMETRICE INTEGRALE PRIME Un sistem de ecuaţii diferenţiale ordinare de ordin superior poate fi transformat intr-un sistem de ecuaţii diferenţiale ordinare de ordin prin introducerea de noi funcţii ca necunoscute Fie sistemul de EDO de ordinul doi cu două funcţii necunoscute:    = = cos , ' y y y x Introducând funcţiile necunoscute auxiliare y , sistemul se poate reduce la forma: Sisteme de ecuaţii diferenţiale ordinare        = = = = cos , , , y y y x y y y y unde avem patru ecuaţii de ordinul întâi şi patru funcţii necunoscute În cele ce urmează, vom considera sisteme de EDO de ordinul cu ecuaţii şi n funcţii necunoscute, sub forma canonică, adică        = = = , , , , n n n n y y y y y y Presupunem că sistemul admite soluţie unică care satisface condiţiile yi yi , ,n Mulţimea        = j = j = j , , , n n n n y C y C y C unde ,C , ,C sunt costante reale, constituie integrala generală sistemului Dacă rezolvăm sistemul în raport cu ,C , ,C atunci integrala generală sistemului are forma        = = = , , , , n n n n F C F C F C Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Funcţiile ,F , , Fn care sunt constante în punctele , se numesc integrale prime ale sistemului Mai general, se numeşte integrală primă sistemului orice funcţie de , care se reduce la constantă dacă se înlocuiesc n , cu funcţii ce constituie soluţie oarecare sistemului Se poate demonstra că orice integrală primă sistemului se poate exprima cu ajutorul funcţiilor ,F , , Sistemul se poate scrie sub forma d , , d , , d , , d x f y y f y y f y y n n n n n = De aceea, în continuare, vom considera sisteme de forma n n n x x P x x P x x , d , d , d cu variabila independentă şi Presupunem că funcţiile , ,n nu se anulează simultan într-un domeniu n Un astfel de sistem se numeşte sistem simetric Definiţia Integralele prime distincte ,F , , ale sistemului se numesc independente dacă matricea i n x F j i , ¶ ¶ are rangul k O metodă comodă pentru rezolvarea sistemului este căutarea unor combinaţii integrabile Pentru găsi combinaţie integrabilă, se poate utiliza următoarea proprietate fracţiilor echivalente: dacă t b a b a b a n = = , atunci, pentru orice ,k , , avem Sisteme de ecuaţii diferenţiale ordinare t k b k a n n n = + + + + Teorema Dacă n m ,m , , sunt funcţii continue pe ÌÂn astfel încât: a + + mnPn pe Ân ; b mn xn pe ÌÂn , atunci D® este integrală primă sistemului Demonstraţie Fie , soluţie a sistemului i n P P x x n i n i , d d Atunci d , + n n n n n x P P P F x şi F unde este constantă ◘ Exemple Să se rezolve sistemul xy z yz y xz x - = d Primele două fracţii formează combinaţie integrabilă Prin integrare, din egalitatea y y x d , obţinem integrala primă C y x Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Din teorema , xy z yxz xyz y y - = + d d Rezultă xy z xy z xyz xy d d , d d = - = şi integrala primă xy În acest fel, am obţinut două integrale prime independente F xy C C y x F z şi sistemul este rezolvat, pentru că , C C x z C x y Să considerăm particulă de masă şi vector de poziţie , într-un câmp de forţe conservativ gradV unde este potenţialul scalar al mişcării Avem ecuaţia diferenţială în m&r& gradV , sau, pe componente, sistemul        ¶ = ¶ ¶ = ¶ ¶ = ¶ , , x V mx x V mx x V mx S && && && Sisteme de ecuaţii diferenţiale ordinare Observăm că x x V x x V x x V m x& &x& x& &x& x& &x& & & ¶ - ¶ ¶ - ¶ ¶ + ¶ , adică t x x t m d d d d + Fie T x& x& x& r& mv energia cinetică particulei de viteză v Atunci V t T t d d d d , de unde deducem că +V const Prin urmare +V este integrală primă a sistemului numită integrala primă energiei APLICAŢII ÎN MECANICĂ, FIZICĂ ŞI INGINERIE Aplicaţia Modele de conflict Modele de tip Lanchester D Comănescu, I Caşu Problema fizică În această aplicaţie vom prezenta modele matematice pentru evoluţia în timp unui conflict între două tabere Notăm dimensiunile celor două tabere, la un moment dat, prin şi pe care le presupunem numere reale În funcţie de tipul de conflict aceste dimensiuni pot reprezenta număr de soldaţi, arsenal militar, etc Modelele matematice prezentate vor fi sisteme diferenţiale de ordinul cu două necunoscute, în care ecuaţiile diferenţiale componente se scriu pe baza unei ecuaţii de bilanţ de următoarea formă: x& PC PO , Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie unde PC reprezintă pierderile în conflict ale taberei PO reprezintă pierderile operaţionale ale taberei datorate bolilor, dezertărilor, etc şi reprezintă reîmprospătarea cu efective taberei x CONCOM simplificat Model matematic Acest model este aplicabil conflictelor convenţionale în care pierderile combatante ale fiecărei tabere sunt presupuse fi direct proporţionale cu efectivul taberei adverse Neglijăm pierderile operaţionale şi convenim că nu există reîmprospătarea celor două tabere Pe baza observaţiilor anterioare, modelul matematic este   = = , , y bx x ay & & unde şi sunt două constante reale strict pozitive Acest sistem diferenţial este liniar şi omogen cu coeficienţi constanţi Ştiind efectivele iniţiale ale celor două tabere obţinem problema Cauchy de evoluţie în timp conflictului        = = = = , , , y y x x y bx x ay & & Soluţie Matricea asociată sistemului este  - - = b a A , având valorile proprii reale şi distincte ab Vectorii proprii corespunzatori sunt , , - b a b a deci soluţia generală ecuaţiei se scrie astfel: Sisteme de ecuaţii diferenţiale ordinare abt abt b k a b k a y t x - = e e , unde ,k sunt două constante arbitrare Aplicând condiţiile iniţiale şi folosind funţiile hiperbolice cosh :Â®Â+ sinh :Â®Â, date de formulele e e , sinh e e cosh z z z z - , soluţia unică problemei Cauchy de mai sus se poate pune sub forma sinh cosh cosh sinh , abt abt a b y x y abt b a x abt = + = + Introducem funcţia :Â ®Â dată de relaţia L bx ay şi arătăm că este lege de conservare: d d × × º ¶ + ¶ ¶ = bx ay ay bx y L x t x L L t t & Interpretare fizică Studiul evoluţiei în timp conflictului poate fi făcut fie prin analiza directă funcţiilor şi fie utilizând legea de conservare L În acest caz vom prefera cea de doua metodă de studiu Din considerente de natură practică ne limităm la primul cadran ÎÂ } Mulţimea { ÎÂ ,L K} este un fragment de hiperbolă în situaţia în care şi semidreaptă ce trece prin origine atunci când Evoluţia în timp conflictului rămâne în mulţimea MK ce trece prin datele iniţiale K bx ay Pentru valori strict pozitive ale constantei se obţine un fragment de Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie hiperbolă situat sub semidreapta definită de iar pentru valori strict negative ale constantei se obţine un fragment de hiperbolă situat deasupra semidreptei definite de K Prin analiza lui şi sau urmărind figura observăm că pentru – situaţia de armistiţiu conflictul nu este câştigat de către niciuna dintre tabere, acestea anulându-şi efectivele simultan Pentru conflictul este câştigat de prima tabără x deoarece efectivele celei de-a doua tabere se anulează mai întâi Pentru conflictul este câştigat de doua tabără y Figura Evoluţia în timp conflictului în modelul CONCOM GUERCOM simplificat Model matematic Acest model este aplicabil conflictelor între două tabere de gherilă în care pierderile combatante ale fiecărei tabere sunt presupuse fi direct proporţionale cu efectivele ambelor tabere Neglijăm pierderile operaţionale şi convenim că nu există reîmprospătarea celor două tabere Modelul matematic este   = = , , y bxy x axy & & Sisteme de ecuaţii diferenţiale ordinare unde şi sunt constante reale strict pozitive Acesta este un sistem de ecuaţii diferenţiale ordinare neliniar, de ordinul I Ştiind efectivele iniţiale ale celor două tabere obţinem problema Cauchy de evoluţie în timp conflictului        = = = = , , , y y x x y bxy x axy & & Soluţie Introducem funcţia :Â ®Â dată de relaţia L bx ay şi arătăm că este lege de conservare: d d × × º ¶ + ¶ ¶ = axy bxy y L x t x L L t t & Utilizăm metoda de reducere sistemului cu ajutorul legilor de conservare De-a lungul soluţiei problemei Cauchy avem integrala primă bx ay bx ay K not - , de unde exprimăm pe în funcţie de obţinând a bx bx ay y Înlocuind pe în prima ecuaţie obţinem problema Cauchy în   , x x x& bx ay bx Ecuaţia diferenţială de mai sus poate fi tratată fie ca ecuaţie cu variabile separabile fie ca ecuaţie Bernoulli vezi cap Împărţind ecuaţia cu obţinem ecuaţia lineară şi omogenă în x u : u& bx ay Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Soluţia generală ecuaţiei omogene asociate este bx ay t uomog , unde este constantă arbitrară soluţie particulară ecuaţiei neomogene este, evident bx ay b u = Deci soluţia generală ecuaţiei lineare în este e bx ay b u bx ay t - = ; rezultă că soluţia generală ecuaţiei Bernoulli de mai sus se scrie astfel c bx ay b bx ay x t - bx ay + = Aplicând condiţia Cauchy, obţinem soluţia problemei Cauchy în x bx ay bx ay t x ay bx x t e - = Înlocuind în expresia lui , obţinem ay bx bx ay t y ay bx y t e - = Expresiile anterioare nu sunt bine definite pentru În acestă situaţie problema Cauchy pentru devine    = = , x x x& bx Aceasta este ecuaţie cu variabile separabile Avem bt x x x x ∫ d , Sisteme de ecuaţii diferenţiale ordinare de unde rezultă soluţia problemei Cauchy: bx t x x t + = Înlocuind în expresia lui , obţinem bx t y y t + = Interpretare fizică În această situaţie, situaţie de armistiţiu, nici una dintre cele două tabere nu câştigă, efectivele ambelor tabere scad, pentru orice moment de timp avem b a y x y t x = iar pentru ®¥ ambele tind la În cazul în care efectivele taberelor combatante scad şi b a x x t®¥ t®¥ ® ® conflictul fiind câştigat de tabăra x În cazul în care efectivele taberelor combatante scad şi a b x y t t ® ®¥ ®¥ , conflictul fiind câştigat de tabăra y VIETNAM simplificat Model matematic Acest model este aplicabil conflictelor între tabără de gherilă x şi tabără convenţională y Pierderile combatante ale taberei de gherilă sunt presupuse a fi direct proporţionale cu efectivele ambelor tabere, iar pierderile taberei convenţionale sunt direct proporţionale doar cu efectivele taberei de gherilă Neglijăm pierderile operaţionale şi convenim că nu există reîmprospătarea celor două tabere Modelul matematic este   = = , , y bx x axy & & unde şi sunt constante reale strict pozitive Acesta este un sistem de ecuaţii diferenţiale ordinare neliniare, de ordinul I Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Figura Evoluţia în timp conflictului în modelul GUERCOM Ştiind efectivele iniţiale ale celor două tabere obţinem problema Cauchy de evoluţie în timp conflictului        = = = = , , , y y x x y bx x axy & & Soluţie Introducem funcţia :Â ®Â dată de relaţia L bx ay şi arătăm că este lege de conservare: d d × º ¶ + ¶ ¶ = axy ay bx y L x t x L L t t & Utilizăm metoda de reducere sistemului cu ajutorul legilor de conservare De-a lungul soluţiei problemei Cauchy avem Sisteme de ecuaţii diferenţiale ordinare bx ay bx ay K not - = , de unde exprimăm pe în funcţie de obţinând b ay ay by x t + = Înlocuind pe în doua ecuaţie, obţinem problema Cauchy în    = = + , y y k y a y& Ecuaţia diferenţială poate fi tratată fie ca ecuaţie cu variabile separabile, fie ca o ecuaţie Riccati Expresia soluţiei depinde de semnul lui K În situaţia în care soluţia este + = y t a y y Înlocuind în expresia lui ăsim + = y t a x x Pentru soluţia problemei Cauchy în este = arctanh tanh y K a t aK a K y , în care am folosit funcţia hiperbolică tanh :Â® , dată de expresia z z z z z - + = e e e e tanh , iar arctanh , ®Â este inversa funcţiei tanh Înlocuind în expresia lui găsim Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie + = arctanh cosh y K a t aK b K x Pentru soluţia problemei Cauchy în este = arctg tg y K a t aK a K y Înlocuind în expresia lui ăsim - = arctg cos y K a t aK b K x Figura Evoluţia în timp conflictului în modelul VIETNAM Interpretare fizică Din considerente de natură practică ne limităm la primul cadran ÎÂ } Mulţimea K { ÎÂ ,L K} M Sisteme de ecuaţii diferenţiale ordinare este un fragment de parabolă Evoluţia în timp conflictului rămâne în mulţimea M ce trece prin datele iniţiale K bx ay Pentru valori strict pozitive ale constantei K se obţine un fragment de parabolă situat deasupra parabolei definite de iar pentru valori strict negative ale constantei se obţine un fragment de parabolă situat sub parabola definită de Prin analiza lui şi sau urmărind figura observăm că pentru situaţia de armistiţiu conflictul nu este câştigat de către niciuna dintre tabere, conflictul durând un timp foarte îndelungat ®¥ Pentru , conflictul este câştigat de tabăra convenţională în timp “foarte mare” ®¥ Pentru K conflictul este câştigat de tabăra de gherilă în timp finit Bătălia de la IWO JIMA Model matematic Bătălia de la Iwo Jima avut loc între SUA şi Imperiul Japonez între februarie şi martie în cadrul campaniei din oceanul Pacific celui de-al Doilea Răzbio Mondial Armata SUA instituit blocadă militară asupra insulei Iwo Jima ocupată de armata japoneză Miza era controlul asupra aeroportului militar de pe insulă Dintr-un total estimat de de militari japonezi prezenţi pe insulă la începutul conflictului peste au fost ucişi iar au fost luaţi prizonieri Insula a fost declarată sigură de către forţele americane în zi conflictului şi toate activităţile de luptă au încetat în zi Interesant este că soldaţi japonezi s-au ascuns în vasta reţea de tuneluri de pe insulă şi s-au predat abia în anul Pentru construcţia modelului matematic se utilizează un model de tip CONCOM în care se neglijează pierderile operaţionale în ambele tabere, armata americană putut să-şi reîmprospăteze efectivele în timp ce armata japoneză nu avut această posibilitate din cauza blocadei Notăm cu efectivele armatei SUA implicate în conflict şi cu y efectivele japoneze Modelul matematic este   = = + , , y bx x ay t & & Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie unde sunt constante strict pozitive iar reprezintă funcţia de reîmprospătare a efectivelor SUA Valorile numerice ale funcţiei de reîmprospătare apar în jurnalul de front ţinut de căpitanul american Morehouse Conform acestui jurnal funcţia are expresia:          £ £ bc a c Mai departe, din aceleaşi notaţii rezultă că > şi deci valoiarea celui de al doilea radical este întotdeauna mai mică decât Astfel, sub primul radical avem totdeauna cantitate negativă şi putem scrie ip unde bc a c p = Luând în considerare formulele lui Euler în particular, eip cos sin , soluţia generală ecuaţiei poate fi scrisă sub forma reală x cos sin cos sin A doua funcţie, poate fi determinată direct din prima ecuaţie sistemului : x k k k x m m x + = && Observând că şi relaţia poate fi adusă la forma x sin + sin iar la forma corespunzătoare x sin sin în care Sisteme de ecuaţii diferenţiale ordinare c p c b a p - = - l , c p c b a p - = - l În al doilea procedeu, vom scrie sistemul sub forma unui sistem de ecuaţii diferenţiale de ordinul introducând două noi funcţii auxiliare necunoscute şi v, , , , v cx cx x v u ax bx x u = = = + = & & & & Conform celor expuse în vom determina valorile proprii ale matricei a sistemului, care satisfac [ det = - l - l - l - l - = c c a b P ceea ce conduce la ecuaţia bipătrată leading to the biquadratic ecuaţia l aceeaşi cu Deducem în acelaşi mod că rădacinile ei sunt pur imaginare, fiind date de Vectorul propriu corespunzător valorii proprii ip este  - - b a p p b a p p i i Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Determinând ceilalţi vectori proprii, corespunzători valorilor proprii ip ip , deducem soluţia generală sistemului sub forma , i i e i i e i i e i i e -i i -i i  - - - - + b  - +  - - - - + a  - =  b a p p b a p p b a p p b a p p b a p p b a p p b a p p b a p p v x u x p t p t unde A+ iB, iD, iar B,C,D sunt constante reale arbitrare Rezultă x Acos Bsin cos sin , b a p A t b a p x cos sin cos sin - - + - = Luând acum arctan D, arctan , obţinem forma soluţiei În final, să observăm că, deoarece problema tratată este problemă de vinraţii, ne puteam aştepta de la bun început la soluţii sub formă trigonometrică Pentru simplificarea calculelor putem deci admite pentru şi expresii de forma sin , sin , = a = a x pt x pt unde B, sunt constante nedeterminate Introducând aceste expresii în sistemul şi simplificând liniile trigonometrice, ajungem la sistemul de ecuaţii algebrice lineare şi omogene Sisteme de ecuaţii diferenţiale ordinare , - = - = Ac p A Bb Soluţia banală defineşte condiţia de echilibru soluţie nebanală va exista doar dacă se anulează determinantul sistemului , = - - c p a b care, dezvoltat, conduce la ecuaţia bipătrată p ; aceasta coincide cu ecuaţia în de rădăcini Deoarece sistemul algebric este omogen, putem determina doar raportul B/ A; calculul corespunzător celor două valori p şi p conduce la şi = cu valorile date precedent Aplicaţia M V Soare , Problema fizică Se consideră coardă verticală întinsă puternic de forţă S Pe coardă sunt fixate, la distanţe egale, trei mase figura Admiţând că, pentru micile deplasări transversale, forţa din coardă nu se modifică apreciabil, se cere să se determine tipurile de vibraţie Model matematic Notăm cu distanţa dintre mase şi cu deplasările transversale ale celor trei mase Ecuaţiile de mişcare se scriu, pentru fiecare dintre particule Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie , , y y a S my y y a S my y y a S my = + = = && && && Figura Coarda întinsă de forţă având fixate trei mase m; ,c ,d cele trei tipuri de vibraţii Soluţie Cu notaţia am S b introducând variabila t bt sistemul capătă forma simplificată d d , d d , d d y y y y y y y y y = t = + t = + t Sisteme de ecuaţii diferenţiale ordinare Metoda matriceală, deşi generală, conduce aici la calcule laborioase De aceea, vom prefera metoda eliminării Observăm că, scăzând ultima ecuaţie din prima, funcţia j satisface ecuaţia diferenţială ordinară lineară şi omogenă de ordinul II, cu coeficienţi constanţi, d d + = t j Ecuaţia caracteristică asociată este l şi deci cos sin cos Mai departe, adunăm prima şi ultima ecuaţie şi deducem d d , d d y y y y y + + t + + t Derivând de două ori ultima ecuaţie şi înlocuind în prima, îl eliminăm pe , obţinând pentru ecuaţia diferenţială de ordinul IV, lineară şi omogenă, cu coeficienţi constanţi d d d d + = t + t y y y Ecuaţia caracteristică corespunzătoare este l cu rădăcinile Deci soluţia generală ecuaţiei este cos cos Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Din cea de doua ecuaţie deducem [ cos cos care, împreună cu determină pe şi [ ] [ cos cos ] cos cos cos , cos + d = d + d = d A A A y A A A y Formulele şi reprezintă soluţia generală sistemului în pentru Aj constante arbitrare Cele trei tipuri de oscilaţii sunt indicate în figura Aplicaţia M V Soare , Problema fizică Să se studieze vibraţiile de translaţie şi rotaţie în planul xOz ale unui bloc de fundaţie aşezat pe un teren elastic Model matematic Ecuaţiile diferenţiale care guvernează fenomenul sunt, în planul Ozx , m&x& kx kxhj Jj&& kj Gh kxh kxhx unde · este momentul de inerţie masic al ansamblului fundaţie-maşină faţă de axa Oy normală pe planul xOz trecând prin centrul de greutate, · greutatea blocului aşezat pe teren elastic, · cota centrului de greutate măsurată faţă de teren, Sisteme de ecuaţii diferenţiale ordinare · deplasarea de translaţie în sensul axei Ox · rotaţia în jurul axei Oy, · kx forţa orizontală produsă de deplasare unitate, · kj momentul în planul xOz produs de rotire unitate figura Soluţie Ecuaţiile şi sunt cuplate; ele pot fi scrise sub forma operatorială d d = + h t m d d = - + kj Gh t kxhx Figura Bloc de fundaţie aşezat pe teren elastic Aplicând operatorul + kx t m d d celei de doua ecuaţii, înmulţind ecuaţia cu kxh şi adunând rezultatele membru cu membru, se elimină deplasarea şi se obţine pentru rotaţia următoarea ecuaţie diferenţială ordinară de ordinul IV, lineară şi omogenă, cu coeficienţi constanţi vezi : Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie d d d d j + = j kj Gh t mk mk mGh Jk t mJ Căutând soluţii de forma exponenţială ert găsim pentru ecuaţia caracteristică = - + + + + j mJ k Gh r m k J k Gh r Mai introducem notaţiile Î + g = + - = j j J mh J J mh k Gh p m k p x x unde px este pulsaţia limită vibraţiei de translaţie în cazul în care nu există rotaţii, iar pj este pulsaţia limită of vibraţiei de rotaţie în absenţa alunecărilor Ecuaţia bipătrată devine gr px pj px pj ; rădăcinile ei sunt obţinute din  + g g = j r px px px p şi sunt toate pur imaginare Vom putea căuta atunci direct soluţia sistemului , sub firma sin , sin , = a j a x pt B pt unde sunt constante ce urmează fi determinate pe baza condiţiilor iniţiale Înlocuind în sistemul şi simplificând liniile trigonometrice, ajungem la sistemul de ecuaţii algebrice lineare şi omogene scris sub forma matriceală =   - - j B C k Gh Jp k mp h x x x Sisteme de ecuaţii diferenţiale ordinare Pentru obţine soluţii nenule, trebuie ca determinanul asociat sistemului să se anuleze; dezvoltându-l, deducem ecuaţia pulsaţiilor p g px pj px pj , care diferă de ecuaţia în doar prin semnul termenului în acest fapt este explicabil prin schimbarea lui în ip Rădăcinile acestei ecuaţii sunt  g g = j p px px px În concluzie, în mişcarea cu două grade de libertate, sistemul maşină-fundaţie poate vibra cu una din pulsaţiile principale sau date de relaţia Raportul amplitudinilor şi ale celor două vibraţii poate fi calculat cu ajutorul primei ecuaţii : p p h p m k h m k k mp k h B C x x x x x x - = - = - = Observaţie Sistemul poate fi rezolvat şi direct, cu metoda matriceală descrisă în scriindu-l mai întâi sub forma unui sistem de ordinul cu ajutorul două funcţii auxiliare şi y , , , j + g y j y = + = j J mh p p p x J hm y hp x y x x x x & & & & Determinantul caracteristic asociat este Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie [ ] g l + g + = + - + g - - l - l - l - j j det p p J mh p p p J hm p hp x x x x x x P Anulându-l, vom găsi valorile proprii ale matricei sistemului care sunt pur imaginare şi coincid cu rădăcinile ecuaţiei în Ţinând seama de forma sistemului, vom putea căuta soluţia sistemului ca şi în cazul metodei precedente, sub forma  + a + a + a + a =  y j D rt C rt B rt A rt y x cos sin cos sin De la acest punct, rezolvarea problemei urmează acelaşi drum ca cel din metoda precedentă Aplicaţia M V Soare , Problema fizică Fundaţia unei maşini de greutate este aşezată pe un mediu elastic figura Aria bazei fundaţiei este iar coeficientul de elasticitate a mediului este ks Pentru evita rezonanţa ce poate să apară în timpul funcţionării, maşina este aşezată pe un cadru rigid, legat de fundaţie prin intermediul unor arcuri a căror constantă elastică echivalentă este Greutatea maşinii şi cadrului este P Se cere să se determine pulsaţiile sistemului fundaţie maşină Date numerice: Q ks m Model matematic Se scriu ecuaţiile diferenţiale ale mişcării măsurând deplasările x şi din poziţia statică de echilibru sistemului; se obţine Sisteme de ecuaţii diferenţiale ordinare m &x& m &x& unde s-a notat ksS Soluţie Sistemul de ecuaţii diferenţiale ordinare lineare de ordinul II , poate fi scris sub forma unui sistem de ordinul I Putem, însă, căuta direct funcţiile necunoscute sub forma exponenţială x ebt ebt ; rezultă , b = b = m C m C Figura Fundaţia unei maşini aşezată pe un cadru rigid Am ajuns la un sistem de ecuaţii algebrice lineare, care are drept necunoscute , C Deoarece acestea nu trebuie să fie simultan nule, rezultă că determinantul asociat sistemului trebuie să se anuleze, adică Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie = - + b D = k k m k , sau + b + m m k k m k k m k Ţinând seama că şi ecuaţia de mai sus devine + b + PQ k Sg Q k k P k g Rădăcinile acestei ecuaţii sunt  - + + + b + PQ k S Q k k P k Q g ksS s , Notăm bi pi Introducând datele numerice, avem , , = = = × × - × × + × × × ×  × + × × × = p de unde , şi Întrucât rădăcinile ecuaţiei caracteristice sunt în acest caz pur imaginare, soluţia sistemului este de forma Sisteme de ecuaţii diferenţiale ordinare sin sin sin sin , a - = = a p t k m p p A k m p x A x t l Aplicaţia Amortizorul dinamic M V Soare , Problema fizică maşină de masă rezemată pe un resort elastic de constantă K, este supusă unei forţe verticale pulsatorii sinwt Deoarece, la anumită viteză de funcţionare maşinii, frecvenţa forţei pulsatorii poate deveni egală cu frecvenţa vibraţiilor proprii ale sistemului existând pericolul de rezonanţă figura a se obişnuieşte să se echipeze instalaţia cu un amortizor dinamic, constând dintr-o masă legată de maşina printr-un resort având constanta elastică k figura Sistemul astfel obţinut are două grade de libertate Se cere să se determine mişcările celor două mase în condiţii iniţiale nule Model matematic Sistemul de ecuaţii diferenţiale ordinare care defineşte deplasările şi la momentul se scrie astfel m&y& M&x& Kx sin wt iar condiţiile iniţiale sunt x = , x& y& Soluţie Cu notaţiile = M K , m k , M k , f M F sistemul devine d d = + x t Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie x t t J - sin d d Figura Rezonanţa sistemului mecanic; amortizorul dinamic Eliminând, ca în aplicaţiile şi funcţia între aceste ecuaţii, găsim t t w  + sin d d În mod asemănător se poate elimina şi funcţia şi rezultă t t t w  + sin d d d d Aşa cum era de aşteptat, operatorul diferenţial aplicat funcţiilor şi este acelaşi, deoarece este linear şi cu coeficienţi constanţi Observând că în ecuaţiile intervin doar derivate de ordin par, putem obţine soluţîi particulare de forma x Asinwt Bsin wt Introducând aceste expresii în deducem Sisteme de ecuaţii diferenţiale ordinare f t N x sinw t N y sinw în care s-a notat N Interpretare fizică Vibraţiile proprii reprezentate de soluţia ecuaţiilor omogene asociate pot fi neglijate, păstrându-se doar vibraţiile forţate reprezentate de soluţia particulară Expresiile arată că masele şi capătă mişcare simplă armonică imediat ce vibraţiile proprii s-au stins Rezistenţa dintre forţa pulsantă şi sistemul se produce atunci când frecvenţa lui şi frecvenţa proprie sistemului sunt egale Luând deci expresiile devin x wt w g = sin , wt w g = sin ; ele demonstrează că amplitudinea lui – care, în mod normal, ar fi infinită se reduce, datorită amortizorului, la valoarea finită f Dacă, în plus, se aleg valorile lui şi ale amortizorului astfel încât w, relaţiile se reduc la x wt g = sin ; ele demonstrează că amortizorul, denumit acordat, anulează în întregime vibraţiile lui M CAPITOLUL STABILITATE Definiţia neformală stabilităţii este proprietatea unui sistem de nu îşi modifica considerabil evoluţia în urma unor mici perturbări ale stării iniţiale Această formulare neriguroasă, care provenit din considerente practice în special din mecanică se traduce în teoria ecuaţiilor diferenţiale prin mai multe concepte, care descriu diverse tipuri de continuitate soluţiei unei ecuaţii diferenţiale ca funcţie de datele iniţiale STABILITATEA SOLUŢIILOR ECUAŢIILOR DIFERENŢIALE Fie , → Ân, [ ∞ Ân care defineşte ecuaţia diferenţială x Presupunem că problema Cauchy definită de ecuaţia şi condiţia iniţială x are soluţie unică ,x ÎD De asemenea, presupunem că ecuaţia are soluţie [ ,¥ ®W Cu notăm norma euclidiană pe Ân Definiţia Soluţia : ,¥ ®W ecuaţiei se numeşte simplu stabilă dacă pentru orice şi orice există > astfel încât pentru orice ÎW cu unica soluţie maximală ,t ecuaţiei care verifică este definită pe [t ,+∞ şi x £e, "tÎ[t ,¥ Sisteme de ecuaţii diferenţiale ordinare b Soluţia [ ,¥ ®Wa ecuaţiei se numeşte uniform stabilă dacă este stabilă şi din proprietatea nu depinde de c Soluţia : ,¥ ®Wa ecuaţiei se numeşte asimptotic stabilă dacă este stabilă şi dacă pentru orice există astfel încât pentru orice ÎW cu unica soluţie maximală , ecuaţiei care verifică x este definită pe [t şi lim ®¥ x t d Soluţia : ,¥ ®W ecuaţiei se numeşte uniform asimptotic stabilă dacă este uniform stabilă şi dacă există astfel încât "t "x cu x unica soluţie maximală , ecuaţiei care verifică x este definită pe [t ∞ şi "e $ > astfel încât "t "x cu x şi orice avem x ,x £e Observaţia Este imediat din definiţia precedentă că orice soluţie a ecuaţiei care este uniform asimptotic stabilă este atât uniform stabilă cât şi asimptotic stabilă De asemenea, orice soluţie uniform sau asimptotic stabilă este simplu stabilă Pe de altă parte, stabilitatea simplă nu implică stabilitatea uniformă, iar stabilitatea uniformă nu implică pe cea uniform asimptotică Observaţia Prin schimbarea de variabilă studiul oricărui tip de stabilitate referitor la soluţia se reduce la studiul aceluiaşi tip de stabilitate referitor la soluţia identic nulă ecuaţiei diferenţiale y’ j’ ; de aceea, se poate admite că t, şi se va studia doar stabilitatea soluţiei identic nule, ecuaţiei Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Considerăm, în continuare, ecuaţiile diferenţiale liniare de forma x unde [ +∞ Mn este funcţie continuă Observaţia Dacă, în general, stabilitatea este proprietate soluţiei şi nu a ecuaţiei, în cazul sistemelor liniare, deoarece prin schimbarea de variabilă , soluţia ecuaţiei se transformă în soluţia identic nulă ecuaţiei , rezultă că dacă soluţia identic nulă ecuaţiei este stabilă respectiv, asimptotic stabilă, uniform stabilă, uniform asimptotic stabilă atunci toate soluţiile ecuaţiei sunt stabile respectiv, asimptotic stabile, uniform stabile, uniform asimptotic stabile Din acest motiv vom vorbi, în continuare, despre stabilitatea ecuaţiei şi vom înţelege stabilitatea soluţiei identic nule sau oricărei soluţii maximale Teorema Următoarele afirmaţii sunt echivalente a Ecuaţia este stabilă b Ecuaţia are un sistem fundamental de soluţii mărginite pe Â+ [ ,+∞ c Toate soluţiile maximale ale ecuaţiei sunt mărginite pe Â+ d Toate matricele fundamentale de soluţii ale ecuaţiei sunt mărginite pe Â+ e Ecuaţia are matrice fundamentală de soluţii mărginită pe Â+ Observaţia Afirmaţii asemănătoare celor din teorema se pot formula şi pentru celelalte tipuri de stabilitate Considerăm, mai departe, ecuaţiile liniare cu coeficienţi constanţi de forma x Ax unde AÎ Mn Matricea se numeşte hurwitziană dacă toate rădăcinile ecuaţiei caracteristice det lIn au partea reală strict negativă Sisteme de ecuaţii diferenţiale ordinare Teorema Dacă ecuaţia este asimptotic stabilă atunci este hurwitziană Dacă matricea este hurwitziană, atunci ecuaţia este uniform asimptotic stabilă Considerăm acum sisteme de ecuaţii diferenţiale ,,perturbate" de forma x Ax unde Mn Â+ → Ân este funcţie vectorială continuă şi local lipschitziană în raport cu al doilea argument, iar pentru orice Denumirea de ecuaţie sistem perturbată ecuaţiei este justificată de considerarea ecuaţiei ca provenind din ecuaţia la care s-a adăugat funcţia perturbatoare , Teorema Fie AÎ Mn Ân, , Â+ W→ Ân o funcţie continuă şi local lipschitziană în raport cu al doilea argument care definesc ecuaţia Dacă există > şi astfel încât exp £M e-wt "t , b " ÎÂ+ W LM "xÎ Ân, Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Propoziţia Forma pătratică este pozitiv definită dacă şi numai dacă există > astfel încât m "x ÎÂn Fie vecinătate deschisă originii în Ân şi Â+ → Ân funcţie continuă pe Â+ şi local lipschitziană pe cu Â+ Definiţia Funcţia Â+ → Â+ se numeşte pozitiv definită pe Â+ W dacă există funcţie Â+ Â+ continuă, crescătoare, cu dacă şi numai dacă r astfel încât V , ÎÂ+ ´W Definiţia Funcţia Â+ → Â+ se numeşte funcţie Liapunov pentru ecuaţia dacă: a este de clasă pe Â+ şi "t Â+, b este pozitiv definită pe Â+ W, c £ ¶ + ¶ ¶ ¶ Σ= t x x V t x t V i n i " Â+ W Teorema Dacă ecuaţia admite funcţie Liapunov atunci soluţia sa identic nulă este simplu stabilă Observaţia Dacă funcţia Liapunov din teorema are proprietăţi suplimentare, atunci se obţin rezultate similare teoremei pentru celelalte tipuri de stabilitate Dacă sistemul este autonom nu depinde in mod explicit de mai precis dacă se consideră sistemul diferenţial x’ unde → Ân, putem căuta funcţia Liapunov independentă de variabila t Sisteme de ecuaţii diferenţiale ordinare Propoziţia Dacă → este continuă pe şi > " x Î atunci există vecinătate originii astfel încât să fie pozitiv definită pe Corolarul Dacă W→ satisface a este de clasă pe şi , b " , c , £ ¶ Σ = x x i n i x V f xÎ W, atunci există vecinătate originii astfel încât să fie funcţie Liapunov pentru ecuaţia autonomă pe Definiţia Fie Rn mulţime deschisă, şi → Ân funcţie de clasă Considerăm sistemul şi definim f’ Sistemul liniar x’ Ax se numeşte liniarizarea sistemului sau prima aproximare acestuia Teorema În condiţiile definiţiei , fie de clasă şi presupunem că f’ este matrice hurwitziană Atunci soluţia banală sistemului este asimptotic stabilă SISTEME DINAMICE AUTONOME Definiţia Fie Rn mulţime deschisă Se numeşte sistem dinamic în o funcţie de clasă : W→ cu proprietăţiile a , " W b " " W Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Definiţia Mulţimea{G : tÎÂ}, se numeşte orbita sistemului dinamic prin se numeşte spaţiul stărilor sau spaţiul fazelor Observaţia Legătura cu ecuaţiile diferenţiale se obţine observând că dacă definim pentru ÎW G f , ¶t = ¶ în ipoteza implicită existenţei derivatei atunci ] ] x t s t t t t s s x G G = ¶ + ¶ ¶ = ¶ ¶ ¶ = = Aşadar, dacă vom avea x’ Definiţia Fie Ân mulţimi deschise Sistemul dinamic din se numeşte topologic echivalent cu sistemul dinamic din dacă există un homeomorfism adică funcţie continuă, bijectivă cu inversa continuă Ân Ân cu h astfel încât orbitele lui sunt aplicate pe orbitele lui păstrând direcţia timpului În această situaţie se spune că şi traiectoriile lui şi sunt topologic echivalente Definiţia Sistemul dinamic autonom parametrizat definit de x’ , xÎ Ân, mÎ Âm se numeşte topologic echivalent în Um Ân cu sistemul dinamic autonom parametrizat definit în Uh de x’ Ân, hÎ Âm dacă a există un homeomorfism în spaţiul parametrilor Âm Âm cu Sisteme de ecuaţii diferenţiale ordinare b "m există hm Ân Ân homeomorfism cu hm Um Vx hm şi care duce orbitele din Um ale sistemului pe orbite din Vx ale sistemului cu h păstrând direcţia timpului Definiţia Fie x’ un sistem dinamic autonom depinzând de parametrul mÎÂ Apariţia, prin variaţia parametrului unor tablouri de fază care nu sunt topologic echivalente se numeşte bifurcaţie COMPORTAMENT PE TERMEN LUNG AL SOLUŢIILOR Definiţia soluţie periodică izolată sistemului se numeşte ciclu limită Prin soluţie periodică se înţelege că există astfel încât pentru orice Â Proprietatea lui de fi soluţie izolată semnifică faptul că există astfel încât pentru orice Ân cu , {x Â} soluţia , nu este periodică Definiţia Considerăm sistemul Se numesc puncte critice sau puncte singulare sau puncte de echilibru punctele xÎ Ân care verifică ecuaţia Definiţia Considerăm sistemul şi fie Ân un punct critic al său a se numeşte punct de atracţie pozitiv invariant dacă există vecinătate a lui astfel încât dacă este soluţie care verifică atunci pentru orice şi lim ®¥ t t Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie b se numeşte punct de atracţie negativ invariant dacă există vecinătate a lui astfel încât dacă este soluţie care verifică atunci pentru orice şi lim ®¥ t t Teorema Fie Ân Ân care defineşte sistemul funcţie de clasă C şi fie Ân un punct critic al lui Dacă este un punct de atracţie pozitiv invariant respectiv, negativ invariant pentru liniarizarea sistemului atunci este un punct de atracţie pozitiv invariant respectiv, negativ invariant pentru sistemul În plus, toate soluţiile care iau într-un arbitrar valoare suficient de apropiată de tind exponenţial la pentru →¥ respectiv, pentru – Teorema Dacă sistemul admite funcţie Liapunov într-o vecinătate a unui punct critic, atunci punctul critic este punct de atracţie pozitiv invariant APLICAŢII ÎN MECANICĂ, FIZICĂ ŞI INGINERIE Aplicaţia Stabilitatea echilibrului pentru pendul A Cernea Considerăm problema oscilaţiilor unui pendul de lungime care se mişcă fără frecare sub influenţa gravitaţiei Să notăm cu spaţiul parcurs de extremitatea liberă a pendulului la momentul şi cu unghiul măsurat în radiani făcut de pendul cu axa verticală la momentul t; avem lx Forţa care acţionează asupra pendulului este F = mg, unde este masa punctului material, iar este acceleraţia gravitaţională Această forţă se descompune după două componente: una având direcţia firului şi care este anulată de rezistenţa firului şi alta având direcţia tangentei la arcul de cerc descris de capătul pendulului Din legea doua lui Newton ecuaţia diferenţială mişcării este mlx” mgsinx, sau, echivalent, Sisteme de ecuaţii diferenţiale ordinare "+ sin x= l g x Dacă se studiază doar oscilaţiile mici, atunci sinx se aproximează prin şi obţinem ecuaţia oscilaţiilor mici ale pendulului "+ x= l g x , care este ecuaţie de ordinul al doilea liniară, cu soluţia generală vezi sin cos t l g t c l g x =c ÎÂ Â Studiem stabilitatea punctelor de echilibru , "t şi =p "t deoarece are semnificaţia de timp, vom lua Arătăm că este simplu stabil, iar este instabil Prin înmulţire cu x’, egalitatea devine x' x"+ x'sin x= g l care conduce, prin integrare pe , t], la cos ' cos x g l x g l Dacă x' pentru şi intr-o vecinătate lui Se obţine, aşadar, că , este punct de echilibru stabil pentru adică este punct de echilibru stabil pentru Acest rezultat este cunoscut ca Principiul Teorema lui Lagrange: Un punct de echilibru al unui sistem conservativ pentru care energia potenţială are un minim local izolat în este stabil Aplicaţia A Cernea Considerăm un circuit electric format dintr-o rezistenţă inductanţă şi un condensator C Notăm cu iR iL iC starea curentului din circuit la momentul unde iR, iL, iC reprezintă curenţii din porţiunile de circuit care conţin rezistenţă inductanţă L şi respectiv condensatorul C Analog, fie vR vL vC starea tensiunilor din circuit la momentul t Din legile lui Kirchhoff se deduce iR iL –iC vR vL vC , iar din legea lui Ohm generalizată g iR vR Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie pentru Din legea lui Faraday se obţine        = = C C L L i t v C v t i L d d d d pentru unde şi sunt inductanţa bobinei şi, respectiv, capacitatea condensatorului C Din aceste relaţii se găseşte că iL şi vC satisfac sistemul de ecuaţii diferenţiale de ordinul întâi        = = , d d , d d L C C L L i t v C v i t i L pentru Pentru simplitate presupunem Notăm iL şi vC Atunci sistemul anterior se scrie sub forma   = = ' , ' , y x x x pentru În plus, presupunând suplimentar că este de clasă derivând prima ecuaţie membru cu membru şi utilizând-o pe cea de doua pentru a-l elimina pe găsim x” g’ x’ pentru Ecuaţia este ecuaţia lui Liénard În cazul în care ecuaţia devine x”+ x’ pentru şi poartă numele de ecuaţia lui Van Der Pol Notăm Sisteme de ecuaţii diferenţiale ordinare - = A şi - = g x b Evident nu este hurwitziană şi, deci, teorema nu se poate aplica În schimb, dacă definim funcţia - - = x y x f , aceasta are liniarizarea în , - - = ' g A Cum această matrice este hurwitziană dacă g’ din teorema rezultă că în acest caz g’ soluţia nulă sistemului de mai sus este asimptotic stabilă Ca consecinţă acestui rezultat deducem că soluţia nulă ecuaţiei Van Der Pol – este asimptotic stabilă Considerăm, în continuare, formă mai generală ecuaţiei Liẻnard x” x’ cu şi xg "x Ecuaţia modelează mişcarea unui oscilator armonic cu frecare Condiţia spune că frecarea este pozitivă, iar condiţia xg indică o forţă de revenire Sistemul canonic asociat este Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie    = = ' ' y x x y Originea este singurul punct critic Fie ∫ x G u d şi y V + , Este uşor de văzut că are un minim în Pe de altă parte ¶ + ¶ ¶ ¶ x x V x y x V În plus, dt d nu se anulează decât în puncte izolate şi este strict descrescătoare pe soluţiile sistemului Aplicând teorema , găsim că originea este punct de atracţie pozitiv invariant Aplicaţia Stabilitatea echilibrului pentru oscilaţii liniare D Comănescu, I Caşu Modelarea matematică fenomenelor fizice care conduc la oscilaţii liniare fost prezentată în aplicaţia “Oscilaţii liniare” de la capitolul Oscilaţii libere Modelul matematic al mişcărilor liniare este: m ·· × × Prin notaţia x · = din ecuaţia diferenţială de ordinul al II-lea se obţine sistemul diferenţial echivalent de ordinul cu necunoscutele : x x x x · ·     × Sisteme de ecuaţii diferenţiale ordinare Se observă că unicul punct de echilibru al sistemului este , Stabilitatea acestui punct de echilibru se face prin metoda directă lui Liapunov utilizând funcţia Liapunov , =w Pentru demonstra că funcţia anterioară este funcţie Liapunov se observă că , este punct de minim absolut strict şi că funcţia este lege de conservare se conservă de-a lungul soluţiilor sistemului Oscilaţii amortizate Modelul matematic al mişcărilor liniare este: m ·· · × × × Prin notaţia x · = din ecuaţia diferenţială de ordinul al II-lea se obţine sistemul diferenţial echivalent de ordinul cu necunoscutele : x x x x m w m · ·    = × ×   Unicul punct de echilibru al sistemului este Stabilitatea se demonstrează utilizând funcţia Liapunov , =w Pentru demonstra că funcţia anterioară este funcţie Liapunov se observă că , este punct de minim absolut strict Derivata în virtutea sistemului funcţiei este: d L L x dt m = × × ×m £ ¶ Urmărind expresiile soluţiilor se vedea aplicaţia “Oscilaţii liniare” constatăm că toate soluţiile tind spre punctul de echilibru când ®¥ ceea ce demonstrează că punctul de echilibru este asimptotic stabil Ecuaţii diferenţiale ordinare cu aplicaţii în mecanică, fizică şi inginerie Oscilaţii liniare neamortizate în câmp gravitaţional Modelul matematic al mişcărilor liniare este: m g ·· × × × Prin notaţia x · = din ecuaţia diferenţială de ordinul al II-lea se obţine sistemul diferenţial echivalent de ordinul cu necunoscutele : x x x · ·     × + Unicul punct de echilibru al sistemului este g w Stabilitatea acestui punct de echilibru se demonstrează cu ajutorul funcţiei g L x w = × Se observă uşor că g L w = şi că punctul de echilibru este punct de minim absolut strict pentru L Derivata în virtutea sistemului funcţiei este: d L L dt x ¶ ¶ = ¶ ¶ L L x x x ¶ × ×w × × = ¶ ¶ , 